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Spatio-Temporal  Analysis 

Archaeology 


What  does  it  mean? 
"ANALYSIS" 


Distinction  and  separation  of  parts  of  a whole  in 
order  to  know  their  principles,  elements,  etc. 

Qualitative  and  quantitative  examination  of 
certain  components  using  specialized  methods 

Study  of  limitations,  characteristics  and  possible 
solutions  to  a problem 


What  does  it  mean? 
SPATIAL  "ANALYSIS" 


Can  we  distinguish  and  separate  parts  of  "space" 
to  get  to  know  its  principles,  elements,  etc..? 

Can  we  examine  qualitatively  and  / or 
quantitatively  the  "spatial  components"  using 
specialized  methods? 


Is  "space"  a problem  whose  possible  solutions  can 
be  formally  addressed? 


What  is  this  thing  we  usually  call 

"space"? 


i is  a relational  structure  called"space"  because  we  move  around  some  things  a 
js  can  move  or  be  moved  around  us 


ve  do  not  move,  or  things  around  us  remain  fixed,  we  hardly  had  an  idea  of  "spc 


"Space"  does  not  exist  as  a physical  thing,  but 
physical  relationships  I have  with  things  around  me 

The  idea  of  DISTANCE 
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in  mathematics,  the  distance  between  two  points  in  Euclidean  space  is 
equivalent  to  the 

segment  length  of  the  straight  line  joining  expressed  numerically. 

In  more  complex  spaces,  as  defined  in  the  non-Euclidean  geometry,  the 
"shortest  path"  between  two  points  is  a straight  segment. 


“place"  (location)  is  NO  INTRINSIC  PROPERTY  TO  OBJECT,  BUT 
FEATURES  OF  THE  OBSERVATIONAL  ACT 


In  general  , in  a physical  system  or  otherwise,  the  term  position  is 
used  to  refer  to  distinct  physical  condition  or  situation  that  the  system 
exhibits  . This  is  common  to  speak  of  the  position  of  the  system  in  a 
diagram  that  illustrates  the  system  state  variables 


“place"  (location)  is  NO  INTRINSIC  PROPERTY  TO  OBJECT,  BUT 
FEATURES  OF  THE  OBSERVATIONAL  ACT 


In  physics  , the  position  of  a particle 
indicates  its  location  in  space  or  space  - 
time.  It  is  represented  by  coordinate 
systems. 

In  classical  mechanics  , the  position  of  a 
particle  in  space  is  represented  as  a 
vector  quantity  with  respect  to  a 
reference  coordinate  system. 

In  general  relativity,  the  position  is  not 
representable  by  a Euclidean  vector,  as 
in  the  space-time  is  curved  in  such  a 
theory  , so  the  position  must 
necessarily  be  represented  by  a set  of 
arbitrary  curvilinear  coordinates  , which 
generally  can  not  be  interpreted  as  the 


In  quantum  mechanics,  the 
representation  of  the  position 
of  a particle  is  even  more 
complex  due  to  non  local 
effects  related  to  the  problem 
of  measurement  of  quantum 

m nr 


What  about  "Time”....? 


...has  it  anything  to  do  with  “space 


“Moment"  (temporal  position)  is  NO  INTRINSIC  PROPERTY  TO 
OBJECT,  BUT 

FEATURES  OF  THE  OBSERVATIONAL  ACT 


A temporal  event  is  an 
expression  of  the  fact 
that  some  kind  of  social 
activity  has  taken  place 
somewhere,  and  that, 
therefore,  there  is  some  f 
characteristic. 


we  also  assume  that  this  location  is  in  a state  s or  the  characteristics 
that  define  that  state  have  changed  or  remain  stationary  (no  change). 


TIME  IS  NO  MORE  THAN  THE  EVIDENCE  OF  CHANGE 
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TIME  IS  NO  MORE  THAN  THE  EVIDENCE  OF  CHANGE 


dimensionality  of  T 

Sequence 

Duration 


Destruction  of  Y 


Construction  of  Y 


Destruction  of  X 


Construction  of  X 


Lifespan  of  X 
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• Isotopic 
events 


Historical  period 


Event  1 (BC)  - Event  4 (BC)  - Duration  of  historical  period 


onsequently,  if  “space"  and  "time"  exists  only  in  our  brain 


What  does  it  mean  "spatio-temporal  analysis" 


searchJD:  ena0108 


Let's  try  to  "analyze"  space 

1 


Can  we  distinguish  and  separate  parts  of 
"space"  to  get  to  know  its  principles, 
elements,  etc..? 


customer 


streets 


parcels 


elevation 


land  usage 


real  world 


The  problem  of  spatial 
decomposition: 

The  definition  of  SPATIAL 
VARIABLES 


Let's  try  to  "analyze"  time 

1 


Can  we  distinguish  and  separate  parts  of 
"time"  to  get  to  know  its  principles, 
elements,  etc..? 

The  problem  of  temporal 
decomposition: 


The  definition  of  temporal 
VARIABLES 


Let's  try  to  "analyze"  time 

1 


Can  we  distinguish  and  separate  parts  of 
"time"  to  get  to  know  its  principles, 
elements,  etc..? 

The  problem  of  temporal 
decomposition: 

The  definition  of  temporal 
VARIABLES 


Let's  try  to  "analyze"  time 

1 


Can  we  distinguish  and  separate  parts  of 
"time"  to  get  to  know  its  principles, 
elements,  etc..? 

The  problem  of  temporal 
decomposition: 

The  definition  of  temporal 
VARIABLES 


Let's  try  to  "analyze"  time 

1 


Can  we  distinguish  and  separate  parts  of 
"time"  to  get  to  know  its  principles, 
elements,  etc..? 

The  problem  of  temporal 
decomposition: 

The  definition  of  temporal 
VARIABLES 


Let's  try  to  "analyze"  time 

1 


Can  we  distinguish  and  separate  parts  of 
"time"  to  get  to  know  its  principles, 
elements,  etc..? 

The  problem  of  temporal 
decomposition: 


The  definition  of  temporal 
VARIABLES 


Let's  try  to  "analyze"  time 

1 


Can  we  distinguish  and  separate  parts  of 
"time"  to  get  to  know  its  principles, 
elements,  etc..? 

The  problem  of  temporal 
decomposition: 

The  definition  of  temporal 
VARIABLES 


Let's  try  to  "analyze"  time 

1 


Can  we  distinguish  and  separate  parts  of 
"time"  to  get  to  know  its  principles, 
elements,  etc..? 


The  problem  of  temporal 
decomposition: 

The  definition  of  temporal 
VARIABLES 


Let's  try  to  "analyze"  time 

1 


Can  we  distinguish  and  separate  parts  of 
"time"  to  get  to  know  its  principles, 
elements,  etc..? 
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The  problem  of  temporal 
decomposition: 


The  definition  of  temporal 
'VARIABLES 


Let's  try  to  "analyze"  time 

1 


Can  we  distinguish  and  separate  parts  of 
"time"  to  get  to  know  its  principles, 
elements,  etc..? 

The  problem  of  temporal 
decomposition: 

The  definition  of  temporal 
VARIABLES 


Let's  try  to  "analyze"  time 

1 


Can  we  distinguish  and  separate  parts  of 
"time"  to  get  to  know  its  principles, 
elements,  etc..? 

The  problem  of  temporal 
decomposition: 

The  definition  of  temporal 
VARIABLES 


THE  FORMATION  AND  DEFORMATION  OF  PREHISTORIC  LAKESIDE 
SETTLEMENT 


FROM  LIFE 
FROM  PAST 


TO  DEATH 
TO  PRESENT 


Depositional 

processes 


Abandonment 

processes 


1 


ARCHEOLOGICAL 

RECORD 


Post  depositional/taphonomic 
processes 


26 


Let's  try  to  "analy 

1 


space-time 


Archaeological 

Space-Time: 

Intersection  of  social 
agents,  social 

activities  and  natural 
processes  in  space 
and  through  time. 

Thousand  of  ways  to 
decompose  space 
into  individualized 
variables: 


Let's  try  to  "analy 

1 


space-time 


Archaeological  Space- 
Time: 

• ABUNDANCE  OF  WELL 
IDENTIFIED 
ARCHAEOLOGICAL 
CATEGORIES:  COUNT 

DATA 

•PRESENCE 
UNCOUNTABLE 
ELEMENTS: 

STRUCTURES 

•ABSENCE:  THE 

IMPORTANCE  OF  EMPTY 
SPAC  E. 


OF 

BUILT 


Let's  try  to  "analyze"  space-time 

1 


Can  we  distinguish  and  separate  parts  of 
"space"  to  get  to  know  its  principles, 


What  is  around  me? 

What  has  been  changing  since  the 
very  beginning? 

Identification  of  observed 
discontinuities. 


Theory  of  Observation:  Physics , 
Chemistry, 

Biology,  Geology. 


Let's  try  to  "analyze"  space*time 

2 


Can  we  examine  qualitatively  and  / or 
quantitatively  the  "spatio-temporal 
components"  using  specialized  methods? 


"Measuring"  space 
implies  defining  and 
quantifying  the 
distance  between 
objects  that  are 
"different". 


Let's  try  to  "analyze"  space*time 

2 


Can  we  examine  qualitatively  and  / or 
quantitatively  the  "spatio-temporal 


comp 


Ain»'>i'hiiii  10  Ex>ri  Bdn*Bf£I00.31  ittO'  c5fl:U  pub  vppkoi] 


1800BC  1600BC  1400BC  1200BC  1000BC  800RC  600BC 


ized  methods? 

"Measuring"  time 
implies  defining  and 
quantifying  the 
distance  between 
two  state  changes 


Calendar  date 


Let's  try  to  "analyze"  space*time 

2 


Can  we  examine  qualitatively  and  / or 
quantitatively  the  "spatial  components"  using 
specialized  methods? 


The  relevance  of  the 
3 dimensions  as  the 
simplest 
representational 
schema.  Below  3 D 
complete 
representation  of 
space  is  impossible. 


Let's  try  to  "analyze"  space*time 

2 


Can  we  examine  qualitatively  and  / or 
quantitatively  the  "spatial  components"  using 


If  we  add  Time  (2 
Dimensions),  the 
minimum  for 
representing 
archaeological 
space-Time  is  5 
Dimensions. 


Representing  Archaeological  Space-Time 


Each  archaeological  event  will  be 
represented  as  a five-dimensional 
vector 


Sj  {x,y,z,t1,t2) 

Formally  speaking,  we  represent 
five  dimensions  as  succesive  four- 
dimensional matrices  (x,y,z,  for 

all  temporal  locations  in  the 
sequences. 


Space,  x 


Let's  try  to  "analyze"  space*time 

3 


Is  "space"  a problem  whose  possible  solutions 
can  be  formally  addressed? 


Problem  Solving  in  Archaeology 


Social  Action Archaeological  Record 

Physical  Properties 
Shape 
Size 

Composition 
Texture 


Locational 

Properties 

Space 

Time 


PRODUCTION 

DISTRIBUTION 

USE 


I 

I 
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Causes 


Effects 


The  Nature  of  the  Archaeological  problem 


Where  activity  X occurred? 


Archaeological  observables  are  not  the  activity,  but  something 
generated  indirectly  during  the  occurrence  of  the  activity 


The  Nature  of  the  Archaeological  problem 


Intuitively:  where  count  data  are  the  more  abundant 


Why  X occurred  HERE  and 
THEN? 


Because  "fire"  happened  there 
And  then 


v 

Why  fire  burned  logs 
there  and  then? 


/ ' f 
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The  research  process 


Explaining  why  social  action 
took  place  there  and  not 
elsewhere 

Simulating  the  results 
of  social  mechanisms 


Measuring  spatial- 
temporal  variability 


Observing 

the 

archaeologic 
al  record 


SICK 


Measuring  space 


Archaeological  Excavation 


Counting 


•Enumeration  of 

archeological  evidences 

within  the  same  variable  or 
unit 


•The  intensity  of  action  (its 
repetitions)  took  place  in  the 
past  within  a spatio-temporal 
unit,  detectable  through  the 
material  consequences  in  the 
archaeological  record 


* • 
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3 4 5 6 
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Frequency  of  archaeological  data: 
an  accumulation  of  some  material  items  on  the  ground 
surface  where  the  action  took  place  or  as  the  intensity  of 
the  action. 


Additional  Difficulties. 

1.  From  Count  data  to  Spatial  Frequencies 


Raw  count  observations  are  not  reliable- 
i/Ve  need  tpo  convert  count  data  into  spatial  frequer 
By  defining  sampling  units  of  equal  spatial  extensio 
and  temporal  duration 


t is  easy  to  stress  extension  uniformity, 
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Figure  P-707 
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When  sampling  areas  are  poligonal  and  irregular 


Calculating  the  centroid  of  a 
polygon  or  the  mean  center  of 
a set  of  points. 


o 


Calculating  the  weighted  mean 
center.  Note  how  it  is  pulled 
toward  the  high  weight  point. 
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(same  example  data  as  for  area  of  polygon) 
n n 

_ 

X=^—,Y=^— 
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Particular  Properties 

Particular  Properties 
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Particular  Properties 

Particular  Properties 
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3 Dimensional  Sampling  Areas 


v count  observations  are  not  reliable- 
need  tpo  convert  count  data  into  spatial  frequer 
defining  sampling  units  of  equal  spatial  extensio 
I temporal  duration 


; not  as  easy  to  stress  duration  uniformity. 

can  use  stratigraphy,  although  it  is  not  always 
liable  as  an  estimate  of  duration 


3 Dimensional  Sampling  Areas 
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Volumetric  calculation.  ROCKWORKS. 


Measuring  Time 
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Stratigraphy 


http://www.harrismatrixcomposer.com/ 


Harris  Matrix  Composer  - HySlraligrjphy.hmo 
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Estratigrafia 

http://www.stratify.org/index.htm 


What  is  Stratify? 


Stratify  1.5 


Home 


Screenshots 


Stratigraphy  is  a method  for  relative  dating  of  archaeological  contexts.  The  Harris  diagran 
visualises  the  relative  chronology  of  a site  with  a diagram  that  shows  the  contexts  and  thei 
stratigraphic  relationships.  The  aim  of  the  program  Stratify  is  to  lay  out  the  Harris  diagrar 
automatically  taking  all  the  available  information  on  chronology  and  groupings  into 
account.  Stratify  stores  the  data  of  the  contexts  and  their  stratigraphic  relationships  in  a 
data  base.  Interactive  and  automatic  checks  are  also  supported  in  order  to  detect 
inconsistent,  missing,  or  contradictory  relationships. 
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data  bases 
Displaying  phases 
Dealing  with  groupings 
Thematic  box  styles 
Term  file  management 
Coordinate  projection 
Interactive  and  automatic 
checks 

Combine  absolute  dates 
and  stratigrgaphic 
relationships  by 
monotone  regression 

Output  formats:  BMP,  JPG,  DXF,  HPGL,  Maplnfo,  SVG,  WMF 
Data  base  import  and  export  formats:  csv,  DBASE 
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Measuring  Time 


Harris  Matrix  nodes  can  explain  only  3 situations 
•Uno  es  after  or  before  another 
•There  is  no  chronological  relationship 
(independence) 

•Both  are  contemporaneous 


representation  of 
structural  entity 


contemporary  with 


earlier  than/later  than 
or  contemporary  with 


earlier  than/later  than 


broadly  contemporary  with 


earlier  than/later  than 
or  broadly 
contemporary  with 
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Measuring  Time 


Temporal  Reasoning  in  an  Arificial  Intelligence 
domain: 

Allen  Algebra  The  ^ow'n®  ^ ^ase  relat'ons  capture  the  possible  relations  between  two  intervals 


Relation 

Illustration 

Interpretation 

X < Y 

X 

X takes  place  before  Y 

Y > X 

Y 

X m Y 

X 

X meets  Y {/  stands  for  inverse) 

Y mi  X 

Y 

lor 

X 

X overlaps  with  Y 

Y oi  X 

Y 

XsY 

X 

X starts  Y 

Y siX 

Y 

XdY 

X 

Y di  X 

Y 

X during  Y 

X f Y 

X 

X finishes  Y 

Y fi  X 

Y 

X = Y 

X 

X is  equal  to  Y 

Y 
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Taphonomy 

A relative  measure  of  Temporal  Duration 


SEQUENCE  OF  PROCESSES 
AFFECTING  PRESERVATION 


SUB-DISCIPLINE 
OF  TAPHONOMY 


ORGANISM  DIES  OR 
SHEDS  BODY  PARTS 

— ORGANIC  (SOFT)  PARTS  DECAY 

SEDIMENTARY  PROCESSES 
INTERACT  WITH  REMAINING  PARTS 

— BURIAL 

CHEMICAL  ALTERATION 
” AND  LITHIFICATION 

— COLLECTION 


NECROLOGY 


BIOSTRATINOMY 


DIAGENESIS 
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3 Dimensional  Sampling  Areas 


Vertical  E*ag  = 1.D0 

v Autodiaw 

qz 

R^ndiy 

L^ 


Fie  Edit  View  Help 


>i'  p-  M ^Stratigraphy  Mott 


Overburden 
Volume  - 57,518,304 
Mass  = 11,503,261 
Aquifer-1 

• Volume  = 25,963,250 
\U^  = 77QB375M 


X~  If  Oriartafen  .'V 

• P -f-  Awes  / 

r@WfldOulinC 
T~  Labels 
Legend: 

Dsla  1 

dS7@  Stratigraphy  Model \ yf 

m □vefburclert  \ 

$ 1^  Ej  Fop  \ 

±J  r MBase 
+J  T @ Sides 

fVoKime  = 57,516_3QA 
Mast -11J503^ 
Aqufei-1 
*J  1^  ®Top 
*J  p 

+i  ^ygSide? 

fVokime  - 25963.250 

Mass  =7,700,975 
Aq^d-I 


n. 


s.aon* 


HI!' 
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MS&Uic 
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Lnked  F-Je  List 


■<53  in1* 


Eiilin 


Overburden  ' od  crd 
DverbiJrcferr_Ra^e  grd 
Aqufet-1_Top  gd 
Aquifct  ’1_B*se.qrd 
Aqutad-1_Tcp.grd 
AqLitafd-l_Base.gid 
Aqufer^Top.cpd 


Layer  Volume  is  independent  of  the  duration  of  formation  process. 


Age  (cal.  AD) 


Calculating  the  time  duration  of  a stratigraphic  unit 

Probability  per  year 

0 0.0 1 0.02  0.03 


Age-Depth  Radiocarbon  Modelling 


14C  An  "absolute"  measure  of  time 


Introduced  in  1946  by  W.  Libby. 


When  the  neutron 
collides™^  a 
ntnoa&n  atom,  a 

nitrogen- 1 A (s(?\rf=n 
protons  'jew-. 
no-iTon* \ a loiim 

luma  into  a 
cerfcon  14 
atom. 


Cosmic  ra  /L  enter  ‘he  earth's 
aimcsphere  ard  coikle  with  cn 
atom,  aealing  an  onorgetio 
neutron. 

* Neuton 


Neurlron  cuplui*? 


L-QrtOn  U 


NrtrOgcr  14 


Proicn 


Plants  ahecrh  carton  dtaxids- 
and  incorporate  carbon-'  4 
Ihnajgh  phctocynthe-sfi 


Animals  and  p&ap1&  ecil 
plarts  ard  take  in 
CCfborvU 


_ FdlowlfiQ  deoih  and 

irst:«i.riol.  wood  end  bone? 


loLe  C- 1 4 oa  1 chants 

to  M- 1 4 by  bftlfa  dea  ':ay 


Beta  decay 


Nitrag^n  1 4 


CarDor  \A 


•'  'oanic  Radiation 
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14C  An  "absolute"  measure  of  time 


Although  mayber  it  is  not  as 
"absolute"  as  it  may  seem. 

It  is  a probabilistic  interval,  and  within 
it  different  times  moments  are 
possible 


UNCALIBRATED  DATE  (BP). 


LABORATORY  CODE  NUMBER  OF  THE  SAMPLE 


"i.  ■-  iKM  I 

ETH-28490  R_Date(32 10,50) 

68.2%  probability 
1522  (68.2%)  1430calBC 
95.4%  probability 
1611  (95.4%)  1406calBC 
Median  1482calBC 


CALIBRATED  DATE  (BC) 


lntCal09  CALIBRATION  CURVE 


68,2%  (la)  PROBABILITY  INTERVAL 


Calibrated  date  (calBC)  9540/,,  (2a)  PROBABILITY  INTERVAL 
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14C  An  "absolute"  measure  of  time 


Measurement  error  . SIMETRIC  (Gaussian,  "normal") 


Normal 


Histograma 


internal  statistical  error  (the  counting 
statistics  of  the  measured  total  counts,  N, 
in  the  series  of  measurements)  and  an 
external  statistical  error  (comparisons  of 
the  error  in  the  mean  of  a series  of  n AMS 
measurements  for  a sample)  have  to  be 
introduced.  The  first  error  is  calculated 
using  the  total  number  of  14C  counts 
measured  for  each  target.  The  second 
one  is  calculated  from  the  reproducibility 
of  multiple  exposures  for  a given  target. 

The  main  source  for  an  assumed  symmetric  (Gaussian)  error  is  produced  in  the  process  of  radiocarbon 
dating  itself,  and  it  is  due  to  the  sample  preparation  in  the  laboratory  and  the  probabilistic  nature  of 
radioactive  decay  measurement.  The  Gaussian  error  is  included  in  the  standard  deviation  associated  to 
the  radiocarbon  date,  as  provided  by  the  laboratories. 


Normal 
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LA  CURVA  DE  CALIBRACION 


UHB 

Univcrsilal  Autonoma 
dc  Barcelona 


La  forma  irregular  de  la 
curva  de  caiibracion 
condiciona  los  resultados 
de  la  caiibracion. 

A un  material  no 
co rres pondera  tan  solo 
una  fecha,  sino  varias 
fee  has  con  dife  rentes 
grades  de  probabiiidad. 
Ademas,  en  algunos 
tramos  la  curva  presenta 
una  forma  plana,  un 
ejemplo  es  la 
"discontinuidad  de 
Hallstatt"  entre  760  y 400 
a.C. 
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i CIENCIES  DE  L’ANTIGUITAT 


COMO  SE  CONSTRUYE  UNA  CURVA  DE 

CALI  BRAC  ION 

Para  corregir  una  datacion  14C  hace  falta  calibrarla. 

Hoy  en  dia  disponemos  de  curvas  de  calibracion  realizadas 
segun  3 metodos: 


Hasta  12000  anos  atras,  datando  anillos 
de  arboles  cuya  edad  ya  esta  conocida 
gracias  a la  dendrocronologfa 


Hasta  30000  anos  atras,  datando  los 
corales  fosiles  a traves  del  metodo  de 
uranio-torio 


Hasta  45000  anos  atras,  datando  las 
varvas  v las  estratificaciones  geologicas 
anuales 


14C  Una  medida  "absoluta”  de  la  Posicion  de  un 
evento. 


!.  Calibration.  NO-SIMETRIC  error  ( non-normal) 


Histograma 


1600,00  1700,00  1800,00  1900,00  2000,00 

nonormal 


The  main  Non-Gaussian  error  is  due  to  the  calibration  curve  and  the  process  of  calibration.  Chemical- 
physical  timescales  (isotopic  degradation)  and  astronomical  timescales  (relative  motion  earth-sun)  are 
not  graded  in  the  same  units;  therefore,  "14C  years"  are  not  necessarily  the  same  as  the  "calendar 
years"  (van  Srtrydonck  et  al.  1999).  This  is  due  because  the  concentration  of  14C  has  not  been  uniform 
all  along  the  time  span  of  the  astronomical  scale  (Aitken  1990,  Bowman  1990) 
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http:  / /ci4.ai’ch.ox.ac.ttk 


ttp : / / b ca  1 . s lief,  ac.uk 


http : / www.  c alpal . d 


C alPal  ( 


Calib  ( 

http : //calib.  qttb.ac.uk  / calib/ 

LA  FIABILIDAD  DE  UNA  DATACION  14C 


Una  datacion  radiocarbonica,  debido  a su  valor  probabilfstico, 
se  caracteriza  por  tener  un  nivel  de  incertidumbre.  Por  lo 
tanto,  hace  falta  reconocer  esta  incertidumbre  a traves  de 

un  rnntrnl  dpi  prrnr  nup  sp  nupdp  intrndurir  pn  la*; 


ERRORES  EN  LA  TOMA  DE  DATOS 


Hay  que  tener  cuidado  con  los  fenomenos  de  bioturbacion  y la 
actividad  humana  que  pueden  desplazar  las  posibles 
muestras. 


Los  meiores  contextos  son: 


» Estratos  en  deposicion  pnmaria  v que  representen  eventos  de  corta 
duracion  (p.ej.  nive!  de  destruccion,  pavimento,  hogar,  etc.) 


Estratos  identificados  por  una  coniunto  ceramlco  sequro 


» Acumulaciones  de  semillas  o de  carbones  en  un  mismo  sitfofp.ej.  en 
una  vasija) 


» Clara  asociacion  de  las  muestras  con  estructuras  mas  grandes  {p.ej. 
carbones  en  un  hogar  o en  un  homo) 


Huesos  en  articulacion.  Los  huesos  largos,  ademas  de  preservar 
mejor  el  colageno,  son  mas  dificiles  de  desplazar. 


, 


ta.  — — — 


LAS  MUESTRAS 


A traves  del  analisis  de  la  cantidad  de  14C  presente  en  los 
restos  organicos  se  pueden  datar  materia les  entre  200  y 
ca.  50000  anos  atras. 


Las  posibles  muestras  se  dividen  en  dos  categorias: 

• Muestras  de  vida  laraa  (carbon  vegetal,  madera) 

• Muestras  de  vida  corta  (semillas.  frutos,  huesos 
humanos  y animales,  dientes,  cuernos) 


Efecto  madera 
antiaua 


Datamos  el  momenta  en  el  que  el  organismo 
interrumpio  el  intercambio  de  14C  con  la 
atmosfera  y no  el  momenta  de  uso  de  la  madera 
que  registramos  a nivel  arqueologico. 


La  fecha  obtenida  es  un  terminus  Dost  auem. 


LAS  MUESTRAS  DE  VIDA  CORTA 


Semillas  v frutos 


Huesos  humanos  v animates  ('colaqeno) 


- >- jf 

Teniendo  un  ciclo  de  vida  mucho  mas  breve 


muestras  de  vida  laraa.  se  caracterizan  por  ser  mucho  mas 
fiables. 


(Tecnicas  de  medicion) 

Actjalmente  existen  dos  tecnicas  para  medir  la  concentracion 
de  14C/12C  de  una  muestra. 


f 

Metodo  de  Libby 


• Metodo  tradicional  para  medir  la 
radiactividad  13  (electrones 
producidos)  de  la  muestra 


R 


/ 


Metodo  AMS 


• Basado  en  la  Espectrometria  de  Masas 
con  Aceleradores  (EMA),  no  detecta  la 
radiacion  emitida  por  los  nucleos,  sino 
que  se  cuentan  directamente  los 
nucleos  de  interes  de  la  muestra. 
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VENTAJAS  DE  UNA  DATACION  AMS 


La  medida  es  mucho  mas  sensible,  se  caracteriza  por  una 
desviacion  estandar  inferior  a la  obtenida  a traves  del 


metodo  tradicional 


Es  mucho  mas  rapida  (ca.  30  m.):  se  mide  directamente 
la  concentracion  de  14C  y no  hace  falta  esperar  a que  los 
atom  os  decaigan 


Las  muestras  son  mil  veces  mas  pequenas  (mg)  que  las 
de  metodo  tradicional,  eso  implica  una  mejoria  en  la 
p re  pa  radon  de  las  muestras  eligiendo  tan  solo  las  zonas 
con  menor  contaminacion 


Oxford  Radiocarbon 
AcceleratajJlit 
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LA  BASE  DE  DATACIONES  14C  DE  CATALUNA 
http://www.telearchaeology.com/cl4 


PORTAOA 


La  base  de  dades  radiocarboniques  de  Catalunya 


BASE  DE  DADES 


BfBUOGRAFIA 


M^TODE  C14 


AGRAIMENTS 


upb  O 

dtBftiwIoim  (te  frttliiiiyl 


Presentacio 


LA  BASE  DE  DADES  RADIOCARBONIQUES  DE 
CATALUNYA  8S  una  recopilacio  de  les  datacions  de 
contexts  arqueoldgics  del  nord-est  de  la  Peninsula 
Iberica.  El  rang  temporal  comenA§a  aproximadament 
50.000  anys  abans  del  present  fins  al  voltant  de  I'any 
500  abans  de  la  nostra  era. 

Aquest  servei  web  posa  a disposicid  de  tothom 
interessat  les  datacions  de  carboni  14,  tal  i com  ban 
estat  publicades,  aba  com  alguna  informacid 
contextual  de  les  mostres  datades  i d'alguns 
materials  arqueoldgics  relevants  que  ban  estat 
associats. 

| veure  mes 


La  base  de  dab 
de  dataciones 
C14  para  la 
prehistoria  de 
Catalunya  esta 
realizada  como 
un  primer  paso 
hacia  un  model 
complejo  de  la 
cadena  de 
inferencia  des  ( 
evento  isotopic 
hasta  evento 
social,  y 
rep  resen  ta  acti 
estado  de 
cuestion. 
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PORTA DA 
BASE  DE  DACES 
BIBLTOGRAFJA 
m£TODE  C 1 4 
EMLLAC05 
AGRA  I HE  NTS 

PFS 


La  base  de  dades  radiocarboniques  de  Catalunya 


Base  de  dades 


Formulari  de  Cerca 


Aquesta  base  de  dades  permet  a Tusuarj  localitzar  mostres  arqueologiques  datades  concretes,  si  es 
conelxen  les  informacions  usuals  (Codi  de  laboratory  Jadment,  etc,),  o be  obtenir  lllstats  de  data  dons 
segons  el  jadment,  la  regio  natural,  i/o  el  peiiode  historic, 

BIBUQGRAFIA 


Mostra 

Q Codi  de  laboradort 
© Material 

Jadment 

Nom 

Terms  municipal 
Com  area 

Q Regio  natural 
Artefacts  opus 

© 

Q Tipus  jadment 
0 Periods 
0 Rank  cronologic 


Ac  an  a I at  a a 

Acanalats  complexos  |£j 

Apendix  bold 
Aurinyacia 


Com  fundona? 

Cap  dels  camps  son 
obligators,  Ef  motor  de 
e#rea  laealitsara  totes 
aquelles  entrades  de  la 
base  de  dades  que 
eolncideiHln  a mb  la 
in  Forma  do  intradulda  per 
I’usuari-  Els  resultats  de  la 
cerca  apareixeran  en 
Format  tabular  mostrant 
valors  da  5 tamps  basin. 
Cad  a liinia  acaba  amb  un 
enllac  ("details")  que 
obra  la  Fitxa  individual  d* 
la  entrada.  Selection  ant  la 
icon  a d'ajuda  a I costat  de 
cad*  camp,  Kumari 
trobara  la  deFinicio  del 
concepter  la  sintaxi 
adopts  da  p*r  a la  baa*  da 
dades  etc,  En  la  recerca 
’avancada1"  es  p-oden  fer 
c*rqu*s  aproximad**- 

o 

Mtouu  u aquHiafia 
de  Cjlilun^i 
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http://www.telearchaeoloqy.com/cl4/ 
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Inaci  Equip  Recerca  Publicacions  Col-laboracions 


Agenda 


Enla^os 


User  login 


You  will  be  redirected  to  the  secure 
CAS  login  page. 


Afegeix  una  datac  a la  DB 


Espanol  I English  I Catala 


Afegeix  tu  Cl 4 dades  a La  base  de  dades  radiocarboniques  de  Catalunya 


Mitjangantaquestformulari  podeu  enviar  les  catacions  radiocarboniques  dels  jaciments  arqueologicS^u^esteu  excavant  a Catalunya.  Les  dades 
enviades  sran  inclosses  a La  base  de  dades  radiocarboniques  de  Catalunya  despres  de  la  revisio  pels  resfJtms^bles  de  la  Base  de  dades. 


Codi  del  laboraton  rtrodu  rt  sense  gons  o espas  (ex.  BetaOOOOO) 

Jaciment:  * 


Municipi:  * 


Afegeix  una  datac  c C 1 4 a la  DB 


Regio  natural: 

i 

STta  utstzat  com  criten  les  8 Reg  crs  Fcrestas  de  Catalunya  (Inventan  Ecc  eg  c i Forestal  de  CataOnya  2008).  Per  veure  la  nomenclatura.  vegeu: 
http  ^/www.te  le  a rc  hae  o lo  gy  .c  o ra  /c  1 4/pf . htn  # re  g io 

Coord  X:* 

i 

A traves  de  este  formulario  podeis  enviar  les  dataciones  radiocarbonicas  de 
los  yacimientos  arqueologicos  que  estais  excavando  a Catalunya.  Los  datos 
enviados  seran  incluidos  a La  base  de  dades  radiocarboniques 
de  Catalunya  despues  de  la  revision  per  responsables  de  la  Base  de  datos. 
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La  cadena  de  la  inferencia 
espacio^emporal  , 


i Eventos  isotopicos  \ i Eventos 
i / i deposicionales 


S\  [ Eventos 

i arqueologicos 
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Eventos  sociales 


FROM  ARCHAEOLOGICAL  INFORMATION  TO  HISTORICAL 


KNOWLEDGE:  a palimpsest  structure 


The  social  action  which 
produced  the  material 
evidence  in  the  archaeological 
record. 


Material  consequences  of  social 
actions  that  happened  in  the  past > 
it  can  be  composed  by  many 
different  particular  depositional 
events , with  different  calendar 
dates  and  different  durations. 


The  particular  action  that  generated 
the  location  of  such  item  at  this 
particular  place  and  moment. 


The  result  of  a radiocarbon  dating: 
the  particular  moment  in  which  a 
living  being  -animal  or  plant- 
ceased  to  interact  with  the 
atmosphere  and  biosphere. 


MUESTRA  DATADA:  Corresponde  a la  materialidad  de  un  suceso  isotopico 
(Coordenadas  cartesianas  locals  (0.0  local) 


CONTEXTO  INMEDIATO:  Corresponde  a la  materialidad  de  un  suceso 
depositacional  y el  suceso  arqueologico  al  que  se  refieren  varios  sucesos 
depositacionales  contemporaneos  relacionados  causalmente.  En  la 
practica  esta  formado  por  los  materiales  arqueologicos  inmediatamente 
vinculados  (espacial  y tafonomicamente)  con  la  muestra  datada.  Unidad 
Minima  de  la  intervencion  arqueologica. 

CONTEXTO  RELACIONADO:  Corresponde  a la  materialidad  de  sucesos 
arqueologicos  causalmente  relacionados  con  el  que  genero  la  muestra, 
tanto  espacial  como  temporalmente.  En  la  practica  esta  formado  por  los 
materiales  arqueologicos  no-inmediatamente  vinculados  (espacial  y 
tafonomicamente)  con  la  muestra  datada,  pero  si  relacionados  de  modo  tal 
que  la  relacion  pueda  expresarse  en  terminos  espaciales  y temporales.  Si 
no  es  posible  la  relacion  espacio-temporal,  el  contexto  secundario  no  sera 
reconstruible.  Asociacion  estratigrafica  de  Unidades  Minimas  de 
intervencion  arqueologica. 


CONTEXTO  GLOBAL  DE  LAS  RELACIONES  ARQUEOLOGICAS  DE  LA 
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MUESTRA  DATADA 
Codigo  Laboratorio 
METODO  DE  DATACION 
Vida  Media  de  la  muestra 
Codigo  Arqueologico  de  Muestra 
Fecha  de  Excavacion 

Fecha  de  Obtencion  y Analisis  de  la  muestra 

Tipo  intervencion  Arqueologica:  excavacion  urgencia,  excavacion  sistematica, 

fondo  de  museo,  prospeccion,  hallazgo  accidental 

Coordenadas  Espaciales  (x,y,z)  en  el  sistema  local  no  georeferenciado. 

Tipo  de  muestra:  Hueso  No  Carbonizado,  hueso  carbonizad,  carbon,  Semilla 
Carbonizada,  Semilla  no  carbonizada  madera-no-carbo, 

Materia  de  la  muestra 

Taxonomia  de  la  muestra 

Porcentaje  Carbon  Marino 

BP  SD  D13C  D13SD  DeltaR  DeltaR  SD 

CALIBRACION:  Curva 

Extremo  inferior  intervalo  Is 
Extremo  superior  intervalo  Is 
Extremo  inferior  intervalo  2s 

Extremo  superior  intervalo  2s 

Referenda  bibliografica  en  la  que  aparece  publicado  el  analisis  cronometrico  de  la 
muestra 

VALORACION  DEL  RESULTADO  CRONOMETRICO  (%).  Calculado  automaticamente 
por  la  Base  de  datos:  100%  = Vida  Media  de  la  Muestra  < 10  anos;  Metodo  de 


Laboratori  d’Arquecrtogia  Quantatativa 


CONTEXTO  PRIMARIO: 

Etiqueta  de  Contexto  primario:  En  el  formato  en  el  que  haya  aparecido  publicada 
(Unidad  Estratigrafica  | Estructura  |) 

LOCALIZACION  DE  LA  MUESTRA  EN  LA  UNIDAD  ARQ.  DE  REFERENCIA:  Ubicacion 
central,  Ubicacion  proxima  al  centro,  Ubicacion  marginal  o periferica 
Tipo  de  Contexto:  hallazgo  aislado,  hallazgo  dentro  de  contenedor 
individualizado,  tumba,  suelo  de  ocupacion,  area  de  actividad,  agujero  de  palo, 
hogar,  silo,  fosa  de  construccion,  matriz  sedimentaria,  acumulacion  no  definida,  ... 
Caracterizacion  del  Contexto:  Cerrado,  Abierto,  Sin  limites  definidos... 

TAFONOMIA:  Primario  (fase  constructiva),  Secundario  (uso  despues  de 
construccion),  Abandono,  Post-depositacional  (amortizacion  despues  del  uso 
primario  y/o  secundario) 

FUNCIONALIDAD  DEL  CONTEXTO  PRIMARIO:  CAMPO  DE  VALOR  UNICO:  Funerario, 
Ritual  No-Funerario,  Residencial,  almacenamiento,  estabulacion  animal, 
produccion  de  utiles,  produccion  de  subsistencias,  defensa,  ... 


ARTE  FACTO -TIPO  RELACIONADO  EN  CONTEXTO  PRIMARIO:  punta  magdaleniense 
IV,  ceramica  cardial,  vaso  campaniforme,  Argar  Forma  6,  punta  palmela,  espada 
de  lengua  de  carpa,  Dragendorff  2^,  Anfora  punica  tipo  Ramon  2.  Sigillatta 
sudgalica,  etc. 


VALORACION  DE  LA  RELACION  CONTEXTUAL  PRIMARIA  (%)].  Valor  calculado  por  la 
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CONTEXTO  SECUNDARIO: 


Etiqueta  de  Contexto  Secundario:  En  el  formato  en  el  que  haya  aparecido  publicada 
(Nivel  Arqueologico  / Fase  de  Ocupacion 
YACIMENTO  ARQUEOLOGICO 

georeferenciacion  se  definira  por  medio  de  cuatro  columnas:  huso,  zona,  x,  y. 
LOCALIZACION  JACIMENT  ARQUEOLOGIC  1:  Termino  Municipal  (Nomenclatura 
Administrativa  Oficial) 

LOCALIZACION  JACIMENT  ARQUEOLOGIC  2:  Region  Natural. 

RELACION  ESTRATIGRAFICA  ENTRE  LOS  CONTEXTOS  PRIMARIOS  QUE  DEFINEN  EL 
CONTEXTO  SECUNDARIO:  DIRECTA,  INDIRECTA.  Si  no  hay  relacion  estratigrafica  de 
contiguidad/contemporaneidad,  el  contexto  primario  no  tendra  vinculacion  a 
contexto  secundario  alguno. 

Tipo  de  Contexto  Secundario:  tumba,  necropolis,  silo,  suelo  de  ocupacion, 
estructura  construida,  area  de  actividad,  poblado,  no  definido 
Funcionalidad  Contexto  Secundario:  CAMPO  DE  VALOR  MULTIPLE  (todas  las 
funciones  identificadas  en  los  distintos  contextos  primarios  que  definen  el 
secundario):  Funerario,  Ritual  No-Funerario,  Residencial,  almacenamiento, 
estabulacion  animal,  produccion  de  utiles,  produccion  de  subsistencias,  defensa,  ... 
ARTEFACTO-TIPO  RELACIONADO  EN  ALGUNO  DE  LOS  CONTEXTOS  PRIMARIOS  QUE 


CONFIGURAN  EL  CONTEXTO  SECUNDARIO:  punta  magdaleniense  IV,  ceramica 
cardial,  vaso  campaniforme,  Argar  Forma  6,  punta  palmela,  espada  de  lengua  de 
carpa,  Dragendorff  2A  Anfora  punica  tipo  Ramon  2.  Sigillatta  sudgalica,  etc. 
VALORACION  DE  LA  RELACION  CONTEXTUAL  SECUNDARIA  (%)].  Valor  calculado  por 
la  base  de  datos  como  resultado  de  la  caracterizacion  del  deposito:  100%= 
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CONTEXTO  HISTORICO: 


PERIODO:  Paleolitico  Superior,  Mesolitico,  Neolitico,  Calcolitico,  Bronce 
Antiguo,  Bronce  Final  lllb,  Hierro  Antiguo  I,  Romano  Altorepublicano,  Baja 
Edad  media  Tardia,... 

CULTURA  ARQUEOLOGICA:  Argar,  vaso  Campaniforme,  megalitismo, 
Sepulcros  de  Fosa,  cardial,  Post-Cardial,  EpiCardial,  Proto  Talayotico,  Romano, 
Mozarabe,.... 

ESTRATEGIAS  ECONOMICAS  REGISTRADAS  DURANTE  EL  PERIODO 
HISTORICO:  CAMPO  DE  VALOR  MULTIPLE:  caza,  recoleccion,  agriculture, 
ganaderia,  metalurgia  (cobre,  bronce,  oro,  hierro),  intercambio,  comercio... 

ESTRATEGIAS  SOCIO-POLlTICAS  REGISTRADAS  DURANTE  EL  PERIODO 
HISTORICO:  CAMPO  DE  VALOR  MULTIPLE:  bandas,  jefaturas,  estado, 
clases_sociales,  despotismo,  imperialismo,  independencia... 
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CONGRESO 


Cronometrias  Para  la  Historia  de 
la  Peninsula  Iberica 


Iber-Crono 

Cronometrias  para  la  historia  de  la  Peninsula  Iberica 


Barcelona,  17-19  de  octubre  2016 


El  congreso  Cronometrias  Para  la  Historia  de  la  Peninsula  Iberica  surge  de  la  necesidad  de 
sistematizar  e integrar  de  manera  exhaustiva  miles  de  dataciones  cronometricas  existentes  y 
sus  contextos  asociados.  Esta  reunion  cientifica  pretende  ofrecer  un  espacio  de  debate  sobre 
temas  teoricos,  metodologicos,  sustantivos  y practicos  acerca  de  la  datacion  de  eventos 
historicos.  Aunque  el  uso  del  metodo  C14  es  mas  frecuente  para  los  contextos  posteriores  a 
50.000  BP,  aqui  queremos  tratar  todos  los  metodos  fisico-quimicos  que  estan  hoy  en  dia  en 
practica,  o en  fase  de  uso  experimental.  Aqui  pensamos  principalmente  en  uranium-torium, 
termoluminiscencia,  geomagnetismo,  dendrocronologia. 


S'esta  connectant  a platform.twitter.com... 


El  C14  y otros  metodos  cronometricos  han  demostrado  su  utilidad  en  la  investigation 
arqueologica  del  pasado  de  la  peninsula  Iberica,  asi  como  en  los  archipelagos  de  la  Baleares, 
Canarias  y Madeira.  Desde  1950  se  han  acumulado  miles  de  dataciones  de  contextos 
arqueologicos.  Sin  embargo,  el  acceso  a las  mismas  es  diffcil,  debido  al  formato  de  su 
presentacion,  a la  inconsistencia  en  la  presentacion  de  los  contextos  asociados  y a la 
dispersion  extrema  de  las  publicaciones.  En  1978  (Almagro-Gorbea  y Fernandez-Miranda)  y 
en  1996  (Lull,  Castro  y Risch)  se  presentaron  esfuerzos  sintetizadores  de  enorme  valor  para 
su  tiempo.  Desde  entonces  se  han  publicado  algunos  trabajos  sinteticos  sobre  periodos 
concretos,  pero  falta  una  vision  general  critica  de  la  potencialidad  de  la  cronometria  para  datar 
la  historia  de  las  sociedades  que  vivieron  en  la  peninsula  Iberica  y archipelagos. 
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Inicio  BASE  DE  DATOS  CURSOS  Inscripcion 


EUBAR:  1567  dataciones  procedentes  de  573  sitios  arqueologicos, 
un  yacimiento  fechado  por  el  14C  cada  916  Km2 


Radiocarbon  determination  (BP) 


Universilat  Auto  noma 
do  Barcelona 


DEPARTAMENTS  de  PREHISTORIA 
i CIENCIES  DE  L’ANTIGUITAT 


r 


y* 


LAv 


Laboratori  dWque^ogia  Quantatativa 


Entry  15 


[ site  | Can  Piteu-Can  Roqueta 

Municipality  Sabadell 
Province/District  Barcelona 
Comarca  Valles  Occidental 
Region  5 

Departament  - 
Canton 

Geographic  Depresion  Prelitoral  (central) 
region 

Country  Espana 

Coord.  X 428122 

Coord.  Y 4598692 

Stratigraphy  - 


[chronology  | 


Laboratory 


KIA-24835 


Year  dating  2004 

Year  BP  (literature)  2755 
Standard  error  30 


Year  BC  recalibrated 
Inferior  boundary  95  978 

Superior  boundary  95  827 
Inferior  boundary  68  Q24 

Superior  boundary  68  842 
Spanish  chronology  Bronce  Final  B 
French  chronology 
Italian  chronology 
German  chronology 


Autonomous  community  Catalunya 


AMS 


[ Function] 


Functional  context  Funerario 

Type  of  site  Necropolis  de  incineracion 


[ Material  | 

Analysed  material  Hueso  humano 
Provenance  Uma  294-34B 


[Context  ] Context  reliability  El  si  □ no 

Settlement  area 

Traces  of  fortification  no 
Elaboration  of  copper  no 
Elaboration  of  bronze  no 
Elaboration  of  iron  no 

Finding  of  copper  no 


Finding  of  bronze  si 

Finding  of  iron  no 

Predominant  domestic  animal 

Evidence  of  agriculture  no  det. 


Necropolis 


Type  of  tomb  Individual 

Funeral  ritual  Incineracior 


Prestige  objects 


Sward  no 

Dagger  no 

Knife  si 

Arrowhead  n0 

Spearhead  no 

Fibula  no 

Pin  no 

Necklace  no 

Earring  no 

Bracelet  no 

Ring  no 

Amber  no 

Ivory  no 

Ax  no 

Horse  no 


C 


Pottery 


Biconical  no 

Carinated  cup  no 

Pottery  with  grooves  si 
Handel  with  vertical  expansion  no 
Decoration  with  “sogueado”  no 
Decoration  with  meanders  no 
Decoration  with  “dientes  de  sierra 
Zig  zag  decoration  si 
Greek  pottery  no 

Phoenician  pottery  no 


SOCIAL  EVENTS  IN  THE  PAST  ARCHAEOLOGICAL  EVENTS  DEPOSITIONAL  EVENTS 
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NEOLITICO  FINAL.  Sepulcros  de  Fosa-Chassey-Lagozza-Cortaillod 
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Spatial  analysis 


The  Nature  of  the  Archaeological  problem 


Where  activity  X occurred? 


Archaeological  observables  are  not  the  activity,  but  something 
generated  indirectly  during  the  occurrence  of  the  activity 


The  Nature  of  the  Archaeological  problem 


Activity  X took  place  there 


Intuitively:  where  count  data  are  the  more  abundant 


An  Intuitive  Model  for  Spatial  Intra-site  Problems 


Where  activity  X occurred? 


probability  an  activity  took  place  somewhere  is  determined  by  the  abundance  and  der 

Of  the  material  consequences  of  that  activity 


An  Intuitive  Model  for  Spatial  Intra-site  Problems 


NELSON'S  Law 


Discard  Equation 
(Schiffer) 


CLARKE'S  effect 


S=  the  number  of  artifacts  normally  in  use 
t=  total  period  of  occupation, 

L = the  use-life  of  the  artifact, 


The  quantity  of  observed  archaeological  remains  is  a function 
of  temporal  duration  of  occupation,  the  size  of  the  involved 
human  group,  and  the  nature  of  the  activity  itself 


Spatial  Prediction  at  the  intra-site  observation  scale 


Is  the  probability  of  an  action 
occurred  at  a particular  place 
the  same  as  the  detection 
probability  of  an  individual  artefact 
linked  to  a particular  activity  at  the 
same  place? 


Spatial  Prediction  at  the  et  ai.  1983 

(see  also  Kintigh  1988,  Shott  1989,....  ) 


the  likelihood  of  an  archaeological  find  within  a given 
survey  area  is  the  product  of  three  probabilities: 

A)  the  probability 
that  an 
archaeologist 
unearths  what  is 
buried  (which  is  a 
function  of  the 
exvacation  area 
size  and  the 
archaeological 
intervention 
layout); 

B)  the  probability  that  an  artifact 
contained  in  the  sampled  area  would  be 
detected  by  the  archaeologist  (which 
among  other  things  depends  on  whether 
or  not  the  depositional  unit  contents  are 
screened,  the  size  and  composition  of 
the  artifact,  and  the  nature  of  sediment); 

C)  the  probability  that  the 

U 1 *^11^ + 


Spatial  Prediction  at  the  intra-site  observation  scale 


the  detection  probability  of  an  individual  artefact  can  be  determined  using  the  following  equation 
(Verhagen  et  al.,  2013) 


d = i - e~AdW 


where 

D = detection  probability; 
e = the  base  of  natural  logarithms  (2.711828); 

A = the  area  of  the  sampling  unit; 
d = the  density  of  artefacts  per  area  unit;  and 

W = the  observation  probability,  i.e.  the  probability  that  an  artefact  will  be  recognized 
such  when  it  is  recovered. 


Spatial  Prediction  at  the  intra-site  observation  scale 


We  can  modify  Verhagen’s  model  of  detection  probability  of  artifacts  adding  Schiffer  discard  equation 


-AdTD 

1-e 


where 

D = Probability  the  action  too  occurred  at  this  place 
e = the  base  of  natural  logarithms  (2.711828); 

A = the  area  of  the  sampling  unit; 
d = the  density  of  artefacts  per  area  unit;  and 

TD  = Number  of  particular  refuse  material  generated  by  that  action  at  that  particular 
place  during  the  occupation  span 


PROCESOS  ESPACIALES 


•Is  observed  clustering  due  mainly  to 
natural  background  variation  in  the 
population  from  which  events  arise? 

•Over  what  spatial  scale  does  any 
clustering  occur? 

•Are  clusters  merely  a result  of  some 
obvious  a priori  heterogeneity  in  the 
region  studied? 

•Are  they  associated  with  proximity 
to  other  specific  features  of 
interest,  such  the  location  of  some 
other  social  action  or  possible  point 
sources  of  important  resources? 

•Are  events  that  aggregate  in  space 
also  clustered  in  time? 


X.5 


Concentrado 


A 


Uniforme 


Aleatorio 


www.blble.ca 


By  Mike  Peters 


"Sorry  about  the  mess,  I invited  a few  of  my 
dog  friends  over  to  watch  the  game  and  things 
got  a little  out  of  hand” 


Conducta  espacialmente  agregada 


Conducta  espacialmente  uniforme 


Conducta  espacialmente  aleatoria 


Aggregated  Spatial  Patterns 


iman  activity  does  not  take  place  everywhere,  but  if  happened  at  some  specific  places  t 
Le  local  conditions  of  places  and  the  social  and  mechanical  needs  of  that  action. 

If  the  theoretical  model  is  correct,  we  should  expect  to  find  aggregated  spatial  patterns,  instead 
of  the  usual  spatial  random  patterns  that  are  used  as  Null-Hypothesis  in  geographical  or 

ecological  models  • 
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Aggregated  Spatial  Patterns 


THE  IMPORTANCE  OF  ZERO  VALUES:  EMPTRY  PLACES. 

Overdispersion:  we  can  say  that  when  points  are  overdispersed,  a point’s  presence  reduces  the 
probability  of  finding  another  nearby 
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Aggregated  Spatial  Patterns 


THE  IMPORTANCE  OF  ZERO  VALUES:  EMPTRY  PLACES. 

Nevertheless,  zero  inflation  may  appear  in  data  sets  with  excess  of  empty  places  caused  by 
false-zero  observations  because  of  sampling  or  observer  errors  in  the  course  of  data  collection. 
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What’s  wrong  with  Poisson? 


It  describes  the  probability  that  a count  unit  appears  at  a 
sampling  (space-time)  unit. 


f(k,  Aj 


k A 


AV 


k\ 


1.  The  probability  of  a success  during  a 
very  small  interval  of  time  At  is  given  by 
a . At. 

2.  The  probability  of  more  than  one  success 
during  such  a small  time  interval  At  is 
negligible. 

3 The  probability  of  success  in  non- 
j^|(yerlapping  intervals  At  are  independent. 


no 


The  Bivariate  Normal  Model 


f{x,y) 

A 


The  Bivariate  Normal  Distribution 


A pair  of  random  variables  X and  I’ have  a bivariate  normal  distribution  iff  their  joint  probability  density  is  given  by 


tf. 


A 2 (1  - p )) 


X - 


^l 


VV  al  ) 


2 2 p (*-  |^)  (y-  n2) 


■ + 


al  a2 


7-M-2 


\ 2 J 


a\ 


2 7T  c?2 


2 

P 


for  -oo  <x<  oo  and  -oo  < y<  oo  , where  > Os  > 0:  and  -1  < p < 1. 


Spatial  Dependence  and  Tobler’s  Law 


Human  dispersal  of  material  elements  fast  never  is  rani 


Everything  is  related  to 
everything  else,  but  near  things 
are  more  related  than  distant 
things. 

-Waldo  Tobler 


112 


Metodologia  de  analisis  de  la  variabilidad  espacial. 


^Como  se  estudia  la  variabilidad  espacial 

observada? 


ESTRUCTURA  CORRELACION 

I 

Determinacion  del  Patron 
Espacial 

Determinacion  del  grado  de 
similitud  entre  localizaciones 


vecinas 


Variation  in  spatial  pattern 
indicative  of  the  INTENSITY  of 
a 

process  at  particular  location 


Kernel  2D  density 


4.  Analisis  espacial  de  poblaciones  teoricas. 

POBLACIONES  TEORICAS 

i.  PROCESO  ALEATORIO  (n=ioo;  n=iooo):  La  ausencia  de  la  relation 

Incremento  en  n provoca: 


- Mayor  irregularidad  en  la  TAS A DE  CAMBIO  = Retfculas  vecinas  sin  relacion. 

- Una  tendencia  a la  homogeneidad  espacial  : 1" 


ii.  PROCESO  ESPACIALMENTE  NORMALIZADO  (n=ioo;  n=iooo):  La 
relation  positiva. 

Incremento  en  n provoca: 

- Reticulas  vecinas  con  relacion. 

- Refuerza  el  caracter  binormal. 
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PROCESO  ESPACIALMENTE  uniforme  (n=325; 
Negativa.  (distribution  Poisson). 


n=i79):  La  relac 


Analisis  espacial  de  poblaciones  teoricas. 

IMPLICACIONES: 

REGISTRO  DE  LA  VARIABILIDAD  ESPACIAL:  ^Coordenados  o Frecuencia? 

Problema:  El  formato  de  datos  condiciona  las  pruebas  a realizar. 

Solution:  Congruencia  en  las  herramientas  analiticas. 

CONSECUENCIA:  La  interpretation  (modelo)  necesita  de  datos  GEOESTADISTICOS 
(frecuencia  espacial),  aunque  lo  ideal  es  partir  desde  COORDENADOS. 


Metodologla  de  analisis  de  la  variabilidad  espacial.  DATOS  EN  COORDENADAS 


Visualization  Patron  Espacial 

KERNEL  DENSITY  ESTIMATION  : 


Figura  35:  Densidad  Kernel  en  2D;  n=100  a la  izquierda,  n=1000  a la  derecha. 


Metodologla  de  analisis  de  la  variabilidad  espacial.  DATOS  EN  COORDENADAS 


Visualizacion  Patron  Espacial 

KERNEL  DENSITY  ESTIMATION 


Figura  35:  Densidad  Kernel  en  2D;  n=100  a la  izquierda,  n=1000  a la  derecha. 


Describing  spatial  variability:  minimum  spanning 

tree 


Metodologia  de  analisis  de  la  variabilidad  espacial.  DATOS  EN  FRECUENCIAS  POR  CUAD 


Visual izaci on  Patron  Espacial 

Diagrama  de  barras  3D  : x / y / frecuencia:spss 


i- 


Where  social  action 
happened? 


parameter 


selection  of 
individual 
observations 


inference  about 
the  population 


point  estimate 


spatial  distribution  of 
productive, 
maintenance  and 
residential  activities 

A 


Archaeological 
intra-site  data 


Where  archaeological 
observables  occur ? 
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FITTING  BIVARIATE  NORMAL 


Test  of  Mardia 


1 "V 

b2,d=  -£[(*,-  IX)'  s~l  (Xi  - fi] 

i= 1 


skeweness  and  kurtosis 
multivariate  test 


Normality  Test 

Equation 

Test  Statistic 

Value 

Prob 

System 

Mardia  Skewness 

14.32 

0. 1538 

Mardia  Kurtosis 

1.5G 

0.1185 

MEASURING  THE  DEGREE  OF  AGGREGATION 


Statistically 
significant 
dispersion 
at  larger 
distances 


rHrHrHrHiMrjrjnjmmnn 


Observed  - Spatial  Pattern 

r Expected  - Random  Spatial  Pattern 

- - Lower  Confidence  Envelop 

— Higher  Confidence  Envelop 


easure  o 
Centrality 
& 


Best  used 

to 

compare 

two 

distributions 

with 

each  other 

May  be  used  for 
polygons  by  first 
obtaining  the 
centroid  of  each 


polygon 


lithics 


Animal 

bones 
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Measure  of 


Identifies  the  mean 
center-point  (X  and 
Y)  of  a group  of 
points. 

Minimizes  sum  of 
squared  distances 
between  itself  and  all 


ts 


Centrality: 

Mean  Center 

Centroid 

Weighted  Mean 
Center 

Center  of  Minimum 


Di 


AT  Y 

x =y — 

tfN 

N v 

y=Y- 


Center  of  minimum 
distance  or  Median 
Center:  Identifies  point  at 
which  distance  to  all  other 
points  is  minimum. 

It  minimizes 

sum  of  distances  between 
itself 

and  all  other  points  Q 


itance 


mm 


in  y_diMD 


Weighted  Mean 
Center:  produced  by 
weighting  each  X and 
Y coordinate  by 
another  variable  (Wi) 

Centroids  derived 
from  polygons  can  be 
weighted  by  any 
characteristic  of  the 
polygon  (as  each 
archaeological 
structural-,, 

Y — ^ 


Mm 


wyi 


IL 


Wi 


Centroid:  The  center  of 
gravity  or  balancing  point 
of  a polygon 

The  center  of  gravity  or 
balancing  point  of  a 
polygon 
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Measures  of  dispersion 


Standard  Distance  Deviation 

Represents  the  standard  deviation  of 
the 

distance  of  each  point  from  the  mean 
center 

Is  the  two  dimensional  equivalent  of 
standard  deviation  for  a single 
variable 

Provides  a single  unit  measure  of  the 
spread  or  dispersion  of  a distribution. 


Standard  Deviational  Ellipse 

It  is  a good  single  measure  of  the 
dispersion  of  the  points  around  the 
mean  center,  but  it  does  not 
capture  any  directional  bias 

The  standard  deviation  ellipse 
gives  dispersion  in  two  dimensions 
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Formulae 

for 

standard 

of 

variable 


Standard  Distance  Deviation 

!• 

Given  by 


1 


Z'jx,-  x.f  +z'JY‘-Y':y 


N 


♦ 

Which  by  Pythagoras: 

l it' ^ 


1 


Or,  with  weights 


^"=iWi(X--  Xcf  +^Wi(Yi-Ycy 


yn  Wi 

i=i 


N 


essentially  the  average  distance  of  points  from  the  center 


We  can  also  calculate  a weighted  standard  distance 
analogous  to  the  weighted  mean  center 


Centrographic  Statistics 


Center  of  Minimum  Distance  or  Median  Center 

Also  called  point  of  minimum  aggregate 
travel 


sdd  = 


i-x.)2+Z^(r,.r,? 


N 


O 


lO 


o 


0 5 10 


i 

X 

Y 

(X  - Xc)2 

(Y-Yc)2 

1 

2 

3 

10.2 

2.0 

2 

4 

7 

1.4 

6.8 

3 

7 

7 

3.2 

6.8 

4 

7 

3 

3.2 

2.0 

5 

6 

2 

0.6 

5.8 

sum 

26 

22 

18.8 

23.2 

Centroid 

5.2 

4.4 

sum  of  sums 

42 

divide  N 

8.4 

sq  rt 

2.90 

Briggs  Henan  University  2010 
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Describing  spatial  variability:  FISRST  ORDER 

ANALYSTS 


Mean  Center  and  Standard  Distance  Deviation 


Mean  Center  and  Standard  Distance  Deviation: 


Sample  size  

Measurement  type 

Start  time  

Unit  


100 

Direct 

06:08:24  PM,  03/22/2007 
Meters 


Variable  

Minimum  

Maximum 

Mean  

Standard  Deviation 
Geometric  Mean  . . . . 
Harmonic  Mean  


X 

1.010000 

6.736000 

4.116740 

0.988640 

3.978663 

3.799889 


3.751000 

9.215000 

6.122780 

1.117231 

6.020372 

5.916382 


Average  Density 
Std  Dist  Dev  . . . 


Circle  Area  Defined 
by  Std  Dist  Dev  . . . . 


3.196230  points  per  sq. 
1.50  m 


7.06  sq  m 


Nearest  neighbor  analysis: 


Directional  mean  and  variance : 


Sample  size : 100 

Start  time : 06:08:24  PM,  03/22/2007 


Mean  angle 

(weighted) 

Circular  variance. 

(weighted) 

Mean  distance 


53.379000 

52.641000 

0.024130 

0.017758 

3.98  m,  13.06  ft,  0.00247  mi 


Intersection  of  mean  angle  (X)  (Y) 

and  mean  distance : 4.203754,  6.124706 

(weighted) : 4.172916,  6.165645 

Triangulated  mean : 3.870396,  5.876943 

(weighted) : 4.116504,  6.122580 


Sample  size 

Measurement  type . . . 
Start  time 


100 

Direct 

06:12:44  PM,  03/22/2007 


Mean  Nearest  Neighbor  Distance  . . : 

Standard  Dev  of  Nearest 

Neighbor  Distance  

Minimum  Distance  

Maximum  Distance  


Based  on  Bounding  Rectangle: 

Area  

Mean  Random  Distance  

Mean  Dispersed  Distance  

Nearest  Neighbor  Index  

Standard  Error  

Test  Statistic  (Z)  

p— value  (one  tail)  

p— value  (two  tail)  


0.18  m 
0.03  m 
6.70  m 


31 .29  sq  m 
0.28  m 
0.60  m 
0.7542 
0.01  m 
-4.7017 
0.0001 
0.0001 


Mean  Nearest 
Neighbor  Distance  (m) 


Expected  Nearest 
Neighbor  Distance  (m) 


Nearest 

Neighbor  Index 


Standard  Deviational  Ellipse: 

Sample  size : 100 

Measurement  type : Direct 

Start  time : 06:08:24  PM,  03/22/2007 

Clockwise  angle  of  Y-axis  rotation..:  39.392497  degrees 

Ratio  of  long  to  short  axis : 2.0837 

SD  along  new  Y axis:  1.91  m 
SD  along  new  X axis : 0 . 92  m 

Y axis  length:  3.82  m 

X axis  length:  1.83  m 

Area  of  ellipse:  5.51  sq  m 

End  time 

Standard  Deviational  Ellipse  (2X) : 


SD  along  new  Y axis:  3.82  i 
SD  along  new  X axis:  1.83  j 


..:  06:08:25  PM,  03/22/2007 


Y axis  length:  7.65  m 

X axis  length:  3.67  m 


Area  of  ellipse:  22.04 


CRIME  III: 

www.icpsr.umich.edu/ crimestat/ 


Point  distribution 


Number  of  points: 

400 

Area  of  convex  hull: 

59,707 

Mean  density: 

6,6994 

Nearest  neighbors: 

Mean  distance: 

0,12241 

Expected  distance: 

0,19316 

Z value: 

-14,017 

p(  random): 

1 ,2275  E-44 

Rvalue: 

0,63365 

Area  estimation 
r Convex  hull 
(*■  Smallest  rectangle! 


Edge  correction  — 
C Off 

(*  Wrap-around 
C Donnelly 


Orientations  and  d stances 


*1 


PAST: 


Describing  data  variability:  convex  hull  and  smallest  rectangle 


Convex  Hull : 

Sample  size  

. . : 100 

Measurement  oype  - - - 

. . : Direct 

Input  units  

. . : Meters 

Start  time  

..:  06:12:44  PM,  03/22/2007 

Points  in  the  convex 

hull:  a 

X 

Y 

4.91700 

4.73100 

4.40000 

3.97300 

3.51400 

3.75100 

1.01000 

4.69000 

3.36000 

7.47000 

5.28500 

9.21500 

6.03000 

9.13000 

6.73600 

7.43100 

CRIME  III: 

www.icpsr.umich.edii/crimestat/ 


# XY  plot 


T3IH 


r LogY 

V 3-pt  average 
r 95%  ellipses 
W Convex  hulls 

V Min.  span  tree 

V Labels 

_xj  Jll  J§| 


<&  Point  distribution 


Number  of  points: 

400 

Area  of  convex  hull: 

35,354 

Mean  density: 

10,423 

Nearest  neighbors: 

Mean  distance: 

0,12241 

Expected  distance: 

0,15453 

Z value: 

-5,012 

p(  random): 

1 . 1 252E-1 5 

Rvalue: 

O', 7906 

| Orientations  and  distances  j 


Area  estimation 
(*■  Convex  hull 
r Smallest  rectangle 

Edge  correction 
C Off 

(*  Wrap-around 
C Donnelly 


PAST 


Testing  spatial  randomness 


A random  distribution  within  its  convex  hull 


Point  distribution 


Number  of  points:  100 

Area  of  convex  hull*  72.66 
Mean  density:  1.3763 


Area  estimation 

(• 

r 

Convex  hull 
Smallest  rectangle 

Nearest  neighbors: 
Mean  distance: 
Expected  distance: 

r 1 1 1 1 1 1 1 1 

^^p(random): 

vtllUb.  — 


0.41045 

0.4262 

■4LZQi09 

0.4795^ 

-ttssITm 


Edge  correction 
r Off 

(*  Wrap-around 
r Donnelly 


% y..v 


— i — 
3.2 


— I — 
4.8 


An  random  distribution  within  its  smallest  rectangle 

D0 


& Point  distribution 


t *•••  •*%% 
• •••  1 • 

Number  of  points: 
Area  of  convex  hull: 
Mean  density: 

100 

76.154 

1.3126 

Area  estimation 
C Convex  hull 
(*  jSmallest  rectangle! 

»*  • ^ • 

Nearest  neighbors: 

. . V 

Mean  distance: 
Expected  distance: 

0 41045 
043642 

Edge  correction 
T Off 

(*  Wrap-around 

• • *••  • • 

c 

p(  random): 

0.2550^ 

C Donnelly 

r • • • ••  • 

k vallJe! 

"TT5405 

Y.5 


Testing  spatial  randomness 


A bivariate  normal  distribution  with  mean=  5,  and 
st.  Dev  = 0.5  within  its  convex  hull 


Point  distribution 


Number  of  points: 

300 

Area  estimation 

Area  of  convex  hull: 

5,7113 

C*1  iConvex  hul( 

Mean  density: 

52,525 

V Smallest  rectangle 

Nearest  neighbors: 

Edge  correction 
P Off 

(*  Wrap-around 

Mean  distance: 

0,064534 

Expected  distance: 

0 ,063983 

"^2, 1 53^. 

C Donnelly 

p(  random): 

0,03239^ 

Rvalue: 

U ,9  3543 

Orientations  and  distances 


*} 


A bivariate  normal  distribution  with  mean  = 5,  and 
st.  Dev=  0.5  within  its  smallest  rectangle 

Point  distribution 


Number  of  points: 

300 

Area  of  convex  hull;: 

5,5501 

Mean  density: 

34,965 

Nearest  neighbors: 

Mean  distance: 

0,064534 

Expected  distance: 

0,054555 

-v,54b/ 

p(  random): 

4.2697E-1 1 

Area  estimation 
P Convex  hull 
(*  Smallest  rectangle! 


Edge  correction  - 
r Off 

(*  Wrap-around 
Donnelly 


Orientations  and  distances 


X | 


Additional  Difficulties 


The  probability  an  activity  took 
place  somewhere  is  determined 
by  the  abundance  and  density 
of  the  material  consequences  of 
that  activity 


From  count  data  to  spatial  frequencies. 
The  problem  of  fragmentation 
The  problem  of  Categorization 
The  problem  of  incomplete  sampling 
The  problem  of  post.depositional  alteration 


Additional  Difficulties. 

4.  The  Problem  of  Incomplete  Surveying 


>v 


Additional  Difficulties. 

4.  The  Problem  of  Incomplete  Surveying 


Sample  Points 


+ 


Sample  Points 

^-4- 


* Unknown  value  "?" 
(to  be  interpolated) 

• "x 


Additional  Difficulties. 

5.  The  Problem  of  Post-Depositional  Alteration 


i 


Pompeii  Premise 


STUDY  OF  INTRA  SITE  SPATIAL  DISTRIBUTION  OF  THE  PAST 

ALTERED  EVIDENCE 


Additional  Difficulties. 

5.  The  Problem  of  Post-Depositional  Alteration 


120  ppm 
1 10  ppm 
ppm 
90  ppm 
BO  ppm 
70  ppfTi 
60  ppm 
50  ppm 
40  ppm 
30  ppm 
20  ppm 
10  ppm 
0 ppm 


3D  Layer  overlaying.  Spatial  Correlation. 


Additional  Difficulties. 

5.  The  Problem  of  Post-Depositional  Alteration 


Ground  surface  correlated  Anisotropy 


Distance 

is  not  3D  interpolation.  It  is  a 3D  view  of  anisotropy 
n different  spatial  gradients 


Influence  of  Ground  Surface  on  Spatial  Distribution 


Additional  Difficulties. 

5.  The  Problem  of  Post-Depositional  Alteration 


Ground  Surface-Abundance  correlation 


3D  Layer  overlaying.  Spatial  Correlation. 


Additional  Difficulties. 

5.  The  Problem  of  Post-Depositional  Alteration 


Orientation  Analysis 


a 


b 


c 


Primary  orientation 
Secondary  orientation 


90  15 


330 


n- 37  cells 


300 
270 

®Hz> 


n=  16  cells 


n- 41  cells 


330 


300 

A.  270  JTt 


Influence  of  Post-Depositional  movement  On  Spatial  Distribution 


Other  methods  for  analyzing  Space  and  Time 


f 


A A A i 

AAA/ 

V ^ A 

AAA/ 

^ A A 

AAA/ 

k A A 

GRADIENT  (1st 
DERIVATIVE) 

Determining  the  rate  of  spatial 
change 

Looking  for  "directionality". 

The  gradient  of  a scalar  field  is  a vector 
field  that  points  in  the  direction  of  the 
greatest  rate  of  increase  of  the  scalar 
field,  and  whose  magnitude  is  that  rate 
of  increase.  In  simple  terms,  the 
variation  in  space  of  any  quantity  can  be 
represented  (e.g.  graphically)  by  a slope. 
The  gradient  represents  the  steepness 
and  direction  of  that  slope. 


Northing  (Feet) 


Gradient  and  directivity 


5 

2 

0 

0 

0 

6 

0 

0 

0 

0 

1 

1 

0 

0 

0 

0 

10 

2 

0 

0 

1 

3 

0 

0 

0 

0 

3 

2 

0 

0 

0 

0 

21 

21 

5 

0 

0 

6 

1 

0 

0 

0 

0 

0 

0 

0 

0 

1 

26 

11 

6 

2 

2 

27 

3 

0 

1 

0 

11 

11 

5 

0 

1 

4 

7 

22 

10 

2 

8 

15 

1 

0 

9 

19 

6 

0 

3 

6 

10 

13 

31 

21 

12 

4 

22 

13 

0 

2 

13 

14 

1 

0 

5 

2 

5 

0 

74 

33 

46 

23 

28 

14 

4 

5 

10 

29 

6 

5 

14 

10 

4 

2 

0 

14 

15 

22 

9 

13 

2 

9 

7 

0 

3 

4 

33 

13 

0 

0 

First  derivative 
surface.  Higher 
areas  are  the 
most  stationary 
and  the  more 
probable  location 
for  slow 
expansion  and 
high  continuity 


SPATIO-TEMPORAL  ANALYSIS  AND  SIMULATION  APPROACHES 


Additional  Difficulties. 

5.  The  Problem  of  Post-Depositional  Alteration 
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The  importance  of  refitting 


A Case  Study.  A Patagonian  Shell 

Midden 


Period  1 Period  2 


Avian  Remains 


A Case  Study.  A Pataagonian  Shell 

Midden 


Sea-wolf  bones  (seals) 
Period  1 

Stationarity  analysis 


Tabla  de  Contingencia  LOBO  MARINO'  EPISODIC  A 


A Case  Study.  A Pataagonian  Shell- 

Midden 


Sea-wolf  bones  (seals) 
Period  2 

Stationarity  analysis 


A Case  Study.  Historical  Archaeology 
A 19th  century  hillfort  in  Central 

Argentina 


Random  pattern 


A Case  Study.  Historical  Archaeology. 
Early  Neolithic  in  Catalonia 


La  Draga  - sector  D - Kernel  density  map  - BOTA,  NAVI 


Random  pattern 


A case  study.  Taphonomic 
Experimentation 


0 100  200  300  400  BOO  0 100  200  300  400  501  0 100  200  300  400  BOO 


ffects  of  scavenging  on  the  frequency  of  bones  (lama  guanicoide).  Controled 
Observations  of  animal  carcasses  during  three  years. 

1st  year:  original  distribution.  Carcass  not  altered  by  scavengers 
2nd  year.  Alteration.  Scavengers  concentrate  some  bones  and  spread  others 

3rd  year.  More  alteration.  Spatial  dispersal 


After:  Mameli,  Barcelo  and  Estevez  2002. 


PREDICTING  THE  ACCUMULATIVE 
CONSEQUENCES  OF  ABANDONMENT 
PROCESSES.  INTRASITE  ANALYSIS  OF 
LAKESIDE  SETTLEMENTS 


ACHINO  Katia  F.,  ANGLE  Micaela,  BARCELO,  Juan 
A. 


laQu 

Laboratori  d'Arqueologia  Quantitativa 


£ GQBIERNO 
)S  DEESPANA 


MINlSTERlO 
□E  ClENCIA 
E INNOVACI6N 


UHB 

Universitat  Autonoma 
de  Barcelona 


SOPRINTENDENZA 
PER  I BENI 
ARCHEOLOGIC1 
DEL  LAZIO 
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VILLAGGIO  DELLE  MACINE(XIX-XVI  B.C.) 


SURFACE  SAMPLING 


MULTIPLE  CAUSAL  MECHANISMS  PRODUCED  THIS 
ARCHAEOLOGICAL  RECORD 


DISECTING 

THE 

PALIMPSEST 


1)  POST-DEPOSITIONAL  ACTIVITIES 

2)  ABANDONMENT 


3)  ORIGINAL  ACTIVITY  DURING  OCCUPATION 
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3)  ORIGINAL  ACTIVITY  DURING  OCCUPATION: 

HOW  TO  PREDICT  THEIR  SPATIAL  LOCATION  FROM  SURFACE  SURVEY  DATA 


PREDICTION  CAN  BE  REALIZED 
THROUGH  EXPLOITATION  OF 
ACCUMULATION 


X 


o-o 

r.FivTRnTn 


Y 


HOW  WE  CAN 
MEASURE  SPATIAL 
ACCUMULATION? 


X 

Y 

MAPPING  THE  SPATIAL 
DISTRIBUTION  OF  THE 
MATERIAL  EVIDENCES 
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INTERPRETING  THE  SPATIAL  DISTRIBUTION 


1 

0 

2 

2 

1 

7 

0 

3 

0 

1 

0 

2 

2 

1 

7 

0 

3 

0 

■ « RANDOMNESS 

HOMOGENEITY  H 


STATISTICAL  REJECTION 
OF  THE  HOMOGENEITY 
HYPOTHESIS 


LOO 


INTERPRETING  THE  SPATIAL  DISTRIBUTION 


POISSON  DISTRIBUTION  AS  NULL  HYPOTHESIS:  ABSENCE  OF  ACCUMULATION 


i 


ABSENCE  OF 
ACCUMULATION 


I 


SPATIAL  HOMOGENEIZATION  jj 
OF  THE  MATERIAL  1 
EVIDENCES  1 

EFFECT  OF  POST-DEPOSITIONAL 
PROCESS  (IN  OUR  CASE-STUDY) 


I 


PRESENCE  OF 
ACCUMULATION 


DEPARTURES  OF  SPATIAL 
HOMOGENEITY 


IT  CAN  BE  EVIDENCE  OF 
ORIGINAL  SPACE  STRUCTURE 
PRESERVATION 


Decreases  ENTROPY 


Increases  ENTROPY 


FAUNAL  REMAINS: 
AREA  2001 


POSTS:  AREA  2001 


PRESENCE/ABSENCE  ANALYSIS 


FAUNAL  REMAINS:  AREA 
2001 


POSTS:  AREA  2001 


NEAREST-NEIGHBOURS  ANALYSIS 


CONVEX-HULL 


Number  of  points:  1 54 

Area:  594,35 

Mean  density:  0,2591 

Nearest  neighbours  analysis  with  wrap-around  edge  correction 
Mean  distance:  0,99508 

Ewpected  distance:  0,98229 

ff:  J ,013 

Z:  0,30909  p (random): 


The  null  hypothesis  of  a random  oattern  (Poisson  orocess) 
can  not  be  rejected  at  ^0.05 


0,75726 


This  is  a consequence  of  exhaustive  survey 
within  the  study  micro-area 


CONVEX-HULL 


Number  of  points:  496 

Area:  810,35 

Mean  density:  0,61206 

Nearest  neighbours  analysis  with  Don  nelly's  edge  correction 
Mean  distance:  0,61556 

Expected  distance:  0,65617 

ff:  0,9351 

Z:  -2,5136  p (random]: 


There  is  statistically  significant  clustering 


0,011951 


This  is  the  result  of  assuming  an  absence  of 
data  at  non  surveyed  areas 
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FAUNAL  REMAINS:  AREA 
2001 


1 .f 

I * »)N> 

\fS 


POSTS:  AREA  2001 


ANALYSIS  OF  SPATIAL  FREQUENCIES  and  not  only  the  density  of  spatial  presence 
location  points 


Kolmogorov-Smirnov  Test 


(^)  = 


Fitted  Distribution 


0,953 


0.I0C 


FAUNA 

Observed  Relative  Frequency 
Expected  Probability 


with  A = 1 
the  spatial 
distribution 
of  our  data 
do  not  fit 
the  Poisson 
distribution 


Kolmogorov-Smirnov  Test 

p-Value 

0 

Fitted  Distribution 


o,:4i 


POSTS 

Observed  Relative  Frequency 
Expected  Probability 


Kolmogorov-Smirnov  Test 

C 0473 

p-Value 

0 

160 


^txan 


Spatial  autocorrelation 

Plot  Mum  Den 


n 


X 


3 p<3.05  criiitai  ewe-' 


1=0 

ABSENCE  OF  SPATIAL 
AUTO  C ORREL  ATION 


Bins  10 


Plot 

* Mcrar'sl 
O Disturce  counts 


THERE  IS  A 
SIGNIFICANT 
DEPARTURE  OF  THE 
HYPOTHESIS  OF 
SPATIAL 
HOMOGENEITY 
MORAN’S  I 


OislarcE 
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Multinormality.  Mardia  Test 


FAUNAL  REMAINS: 
AREA  2001 


Mardia  tests 

Parameter  Value 

Skewness:  47,3 

Skewness,  small 
sample  corrected: 

Kurtosis:  73,33 

Doornik  and  Mansen  omnibus 

FpC  1.521E04 

p (normal): 


POSTS:  AREA  2001 


n 

b2,d  = \ y KXi  - /!)'  s-1  (Xi  - n)] 

n L^l 


1=1 


Statistic 

3337 

8373 

175,4 


df 

10 

10 


P (normal) 

0 

0 

0 


I 


Mardia  tests 

Parameter 
Skewness: 

Skewness,  small 
sample  corrected: 

Kurtosis: 

Doornik  and  Hansen  omnibus 

Fri  5357 


SIGNIFICANT  DEPARTURE  OF  THE  BI- 
NORMALITY’S CONDITION: 
THERE  IS  NOT  A SINGLE  INTENTIONAL 
ACCUMULATION 


Value 

Statistic 

df 

£■  (normal} 

1,631 

303,4 

10 

2,9  E- 59 

304,7 

10 

l,563E-59 

13,63 

-3,574 

7,074  E-05 

FAUNAL  REMAINS: 
AREA  2001 


Grldding  (spatial  interpolation) 


Settings  Semivariogram  Plat  Numbers 


Crass  validation 

Model: 


GaLssiar 


Nugget:  1,2047 

Range: 

Scale: 


SO,3610 


145566 


Optimize 


SSerror:  10,771 

Bins:  10 


Plot 

5emivariogram 
Distance  counts 


Distance 


1) Small  nugge  effect 

2)  Function  increasing  at  the  asymptote:  when  the 
distance  h increase,  also  increase  the  variance.  Two 
accumulations  in  two  far  from  areas,  in  the  two  extremes 
of  the  investigated  area. 

3) Medium  level  of  spatial  irregularity  due  to 
fragmentation 

4) Classical  model  of  accumulation:  majority  of  finds  in 
two  areas  where  the  accumulation  activity  took  place  and 
in  the  surrounding. 


SPATIAL  PREDICTIONS 


Kriging 


Differential  accumulations  of  faunal  remains  and 
a significant  empty  area  between  both 
accumulations 
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Griddfng  (spatial  interpolation) 


Settings  Semivanogram  Plot  Numbers  Cross  validation 

Model: 

GaLssian  v 


Distance 


2) The  population  is  distributed  according  to  a 
more  standard  function. 

3) The  higher  values  are  concentrated  in  the  first 
portion  of  the  function  (from  0 to  24)  and  then 
they  decrease. 

4) The  variance  decrease  in  correspondence  with 
the  increase  of  the  distance  between  points. 

5) The  function  shows  a linearity  (posts 
alignments):  spatial  connectivity  due  to  the  partial 
preservation  of  ancient  built  structures 


SPATIAL  PREDICTIONS 


FAUNAL  REMAINS 
AREA  2009 


Posts:  AREA  2009 


Multinormality.  Mardia  Test 


n 

b2,d  = \ y [(*«  - /O'  5"1  (Xi  - n )] 

n L—k 

i= 1 


Mardia  tests 


Parameter  Value 

Statistic 

df 

^{normal} 

Skewness:  13,35 

Skewness,  small 

2966 

10 

0 

sample  corrected: 

2976 

10 

0 

Kurtosis:  31,19 

53,96 

0 

Mardia  tests 


Parameter  Value 

Statistic 

df 

£■  {normal} 

Skewness:  5,332 

Skewness,  small 

1279 

10 

1.12E-268 

sample  corrected: 

1284 

10 

1.249E-269 

Kurtosis:  18,7 

12,19 

0 

Doom i It  and  Hansen  omnibus 

3606 

^{normal}: 


Y 


Doornik  and  Hansen  omnibus 

2682 

^{normal}: 


SIGNIFICANT  DEPARTURE  OF  THE  BI- 
NORMALITY’S CONDITION: 
THERE  IS  NOT  A SINGLE  INTENTIONAL 
ACCUMULATION 


FAUNAL  REMAINS: 
AREA  2009 


^ Gridding  (spatial  interpolation) 


Settings  Sermivariogram  Plot  Numbers  Cross  validation 


Model: 

Cubic 


Nugget:  OrQDOO 


Range:  42r62fX) 


Scale:  7,7339 


Optimize 
SSeiror:  47,136 

Bins:  10 

Plot 

° Semi  venogram 
Distance  counts 


1)  Increasing  monotonic  function:  when  the 
distance  h increase,  also  the  variance  increase. 

2)  The  function  start  from  the  0-0  but  the  points 
are  attested  from  4.8. 

3)  An  area  of  faunal  remains7  accumulation  is 
recognized  while  the  surrounding  shows  high 
rate  of  empty  areas:  classical  model  of 
accumulation. 


SPATIAL  PREDICTIONS 


Kriging 


Differential  accumulations  of  faunal  remains  and  a 
significant  empty  area  between  both  accumulations 
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Posts:  AREA  2009 


Gridding  (spatial  interpolation) 


Settings  Semivariogram  Plot  Numbers  Cross  validation 


Model: 


GaLssian  v 

Nugget: 

Or21 1 2 

Range: 

6S..11SO 

Scale: 

OrS149 

Optimize 

SSeiror: 

0{ 075526 

Sins: 

10 

Plot 

Semivariogram 
O Distance  counts 


Distance 


1) Small  nugge  effect  (from  0 to  0,20). 

2) The  population  is  distributed  according  to  a more 
standard  function. 

3) The  higher  values  are  concentrated  in  the  first 
portion  of  the  function  (from  0 to  15)  and  then  they 
decrease. 

4) The  variance  decrease  in  correspondence  with  the 
increase  of  the  distance  between  points. 

5) The  function  shows  a linearity  (posts  alignments) 

6) One  outlier  is  also  attested. 


SPATIAL  PREDICTIONS 


Kriging 


Differential  accumulations  of  posts  and  a significant 
empty  area  between  both  accumulations 


Micro  / Macro 


Micro:  Supuesto  Inicial:  Espacio  Continuo 
Macro:  Supuesto  Inicial:  Espacio  Discreto 


Other  case  studies:  Urban  processes  in  Early  Middle  Ages 


100  spatio- 

temporal 

units 

We  studied  the 
diachronic 
changes  in 
density  and 
dispersion 
patterns  around 
an  Early 
Medieval  city. 
Our  goal  is  to 
understand  the 
urban  influence 
in  the 

Islamisation 

process. 


POSSIBLE  SOLUTIONS 

Taking  spatial  heterogeneity  into  account 


TRADITIONAL 
MODELS 

Von  Thunen,  Weber, 
Christaller, 
Hagerstrand 


PROCESSES  ARE  NOT 
BY  SPATIAL  LOCATION 


( ^ 

SPATIAL 

HETEROGENEITY 

MODEL 


THE  PROBABILITY  OF  DIFFERENT  ACTIONS 
VARIED  THROUGH  TIME  AND  SPACE 


Spatial  Heterogeneity  & Anisotropic  Distances 


Spatial 

Heterogeneity 


Directional 
Dependence  of 
Movement 


Use  of  Non- 
Euclidean 
Distances 


Spatial  Heterogeneity  & Anisotropic  Distances 


Spatial 

Heterogeneity 


Directional 
Dependence  of 
Movement 


\ Use  of  Non- 
) Euclidean 

/ Distances 


The  Problem 


Llegenda 

0 Tortosa 
• assentaments 


Extract_dist3 

<VALUE> 

■ 0 - 733 


2.201  - 2.935 


2.935  - 3.669 


3.669  - 4.403 


4.403-5.137 


0 2.500  5.000  10.000  Meters 

1 i i i I i i i I 


DENSITY  (x2TEST) 


DISPERSION  (LCN) 


Why  activities  occurred  where  they  did? 


Why  activities  occurred  where  they  did? 

▼ 

Because  other  people  did 
what  they  did  nearby 


Why  activities  occurred  where  they  did? 


t 


Because  other  people  did 
what  they  did  nearby 


I’M  NOT  DOING  WHAT  MY  NEIGHBOUR  DOES 

HOWEVER,  MY  RATIONAL  AND  IRRATIONAL  DECISSIONS 
ARE  AFFECTED  BY  THE  WAY  MY  NEIGHBOURS 
TAKE  RATIONAL  AND  IRRATIONAL  DECISSIONS 


Why  the  fire  occurred  here? 


Spatial  Correlation 


W/M/ 


QrtPuiMap 


Iq..  iJj-.l-iii 
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Correlating  Differentiated  Spatial 

Distributions 
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Spatially  Constrained  Correspondence 

Analysis 


IANA5HUAIA 


1 ° I ^ li 


File  Edit  Transform  Plot  Statistics  Multivar  Model  Diversity  Time  Geomet  Strat  Cladistics  Script 
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11 


12 


13 
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IE- 


16 
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IS 


19 
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21 
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23 


24 


25 


26 


Edit  made  Edit  labels 


r Square  mode 


□ 
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3669, 
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3669, 

4505, 
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4505, 

4505, 

4505, 

-505. 

4505, 

4505, 

4505, 

-EOE. 

4505, 

4505, 

5341. 


5-341 
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6590,54 

9457,63 

10324,71 

11191,79 

12056,66 

12925,96 

13793,04 

14660,13 

15527,21 

16394,29 

17261,36 

16126,46 

6590,54 

9457,63 

10324,71 

11191,79 

12056,66 

12925,96 

13793,04 

14660,13 

15527,21 

16394,29 

17261,36 

16126,46 

6590,54 

9457,63 
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Spatially  Constrained  Correspondence 

Analysis 


Axe  2:  35  % inertia 

0.3 -i 
0.6- 


* Pottery  1 


1:  56%  inertia 


■0.2- 


— I 1 1 — 

-1.3  -1.5  -L2 


“i r — 

-0.9  -0.6 


1 — 

-0.3 


0.3  0.6 


-0.2- 
-04- 
-0.6- 
-0.3- 

1 Animal  Bone 


Axis  1 Axis  2 Axis  3 


•Lithic 


1 -0.293896  -2.201040.76671 

3 -1.602281.67154  2.4952 

4 -1.546960.931213  -6.47682 

6 -1.08506-0.0227083  - 

0.709457 

8 -0.906616  -0.535651  - 

0.710364 

9 -0.0466913  -1.77221 

0.625382 

13  1.43653  0.800765 
0.222589 

14  1.43653  0.800765 
0.222589 

15  0.571319  -0.700135 

0.27206 

16  -0.124104  0.22262  - 

0.825263 

17  0.200513  -1.34338 

n zLSzLHSS 


Spatially  Constrained  Correspondence 

Analysis 


X 

Y 

Axis  1 
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3 * 
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In  red  locations  where  Pottery  1 and  2 is  abundant 
In  blue  locations  where  Pottery  1 and  2 is  scarce 


Space-Temporal  Analysis. 
An  example 
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SPACE-TEMPORAL  ANALYSIS  OF 
RADIOCARBON  EVIDENCE 
AND  ASSOCIATED  ARCHAEOLOGICAL 

RECORD: 

FROM  DANUBE  TO  EBRO  RIVERS 
AND  FROM  BRONZE  TO  IRON  AGES 


Giacomo  CAPUZZO,  J.A.  BARCELO 


Department  of  Prehistory,  Autonomous  University  of  Barcelona 


Expansive  phenomena 


if  Asiatics 

mi'i- 


system  such  that  every  point  of  the 
?■  e?/?-  r space  has  a distinctive  behavior  through 

^ - w-  Jit'  < .-/nav .*  . r r 


WM  WW  - 


w&>  " 


3 in  volume  resulting  from  an  increase  in 


phy 

n which  a gradient  of  a scalar  field  can  be 


with  the  notion  of  directivity 
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The  historical  problem 


"Urnfield  culture"  (1300-650 
BC).  A process  of  Cultural 
Standardization 


o 

o 

0 


Around  1100  BC  Europe  was  characterized  by  a 
number  of  quite  distinct  regional  traditions  rooted  in  or 
influenced  in  the  Urnfield  tradition  of  shared  burial 
rituals  and  religion. 

The  diffusion  of  cremation  burials  at  the  end  of  the 
Bronze  Age  has  been  traditionally  seen  as  an  expansion 
from  the  Danube-Carpathian  area  toward  northern  and 
western  territories. 

The  possibility  of  expansive  processes  has  been 
traditionally  detected  also  for  other  variables,  like 
specific  pottery  typologies. 


How  can  we  quantify  such  phenomena? 
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Expansive  phenomena: 
European  Bronze  Age  as  a case  study 

> Population  Substitution:  Invasion,  migration,  colonization 
(Wave  of  Advance) 


> Adoption  of  a "foreign"  idea  (Cultural  Transmission, 
acculturation).  "Knowledge"  or  "behavior"  travel  even  in  the  case 
people  do  not  move. 


> Social  Influence  Theory 

(adoption  of  innovation) 


© Original  Artist 
RetoroductjonTighOs’obtajnable  from 
www.CartoonStock.com 
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Estamos  en  el  ano  50  antes  de  Jesucr 
Galia  esta  ocupada  por  los  romanos 
Una  aldea  poblada  por  irreductibles  gale 
davia  y siempre  al  invasor.  Y la  vida  no 
las  guarniciones  de  legionarios 
cidos  campamentos  de  Babaorum. 
num  y Petibdnum... 
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FROM  EAST  TO  WEST  AND  THE  OTHER 

WAY  AROUND 


Urnfield  Culture  1300  -7*  0 B 
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EUBAR:  ^C-dated  archaeological  contexts 


It  includes  more  than  1700  14C 
dates  with  information  about  a 
large  variety  of  associated 
archaeological  contexts  from 
the  Danube  to  the  Ebro  rivers 
and  from  the  Bronze  to  the  Iron 
Ages. 


Types  of  contexts 


■ Residential 

■ Funerary/cultual 

■ Mine 

■ Smithing  site 

■ Infrastructure 

■ n.a. 


Years  BC 


The  traditional  conventional  chronology 


N of  the  Alps  S France 


N-E  Iberian 

N Italy  n . 
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COMPARACI  N DE  LAS  DATACIONES  DE  DIFERENTES  CONTEXTOS 
ARQUEOL  GICOS  DE  CATALU  A (Se  han  eliminado  las  dataciones  con 

desviaciones  t picas  superiores  a 80  adGOS)  (Fuente:  BARCEL  J.  a.  2007-08,  La  seq  tncia  crono-cultural  de  la 
prehistoria  catalana.  Analisi  estad  stica  de  les  datacions  radiomttriques.  Cypsela  17:65-88.) 
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Laboratori  d’Arqueotogia 


Sequential  model  for  archaeological  contexts  located  in 
Northern  Italy  (Amode|=98.4) 


OxCal  v4  2.3  Bronk  Ramsey  (2013);  5 IntCal  13  atmospheric  curve  (Reimer  et  al  2013) 
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Radiocarbon  dates  (1  a)  from  reliable  archaeological 
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CRONOLOG  A DE  LA  EDAD  DEL  BRONCE  EN  EL  SUR  DE  FRANCIA 
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Transition  boundaries  of  the  contiguous  model  for  archaeological 
contexts  located  in  Southern  France  (Amodc,= 145.7) 
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RESULTS  OF  THE  BAYESIAN  MODELING  FOR  NORTHERN  ITALY 
AND  SOUTHERN  FRANCE.  Only  the  analyzed  sub  phases  have  been 
represented.  The  conventional  chronology  is  shown  above  the  x-axis. 
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INFERRING  POPULATION 
GROWTH  FROM  14C  DATA 
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THE  STUDY  OF  POPULATION  TRENDS  IN  THE 
BRONZE  AGE  AND  IN  THE  IRON  AGE  TRANSITION 


Kristiansen  (1998b)  suggested  that  the  Urnfield  period  was  one  on  the  most 
densely  populated  phase  in  the  Prehistory  characterized  by  an  enhance  of 
population  size  in  a large  part  of  Europe,  mainly  in  the  period  1100/1000  BC. 


For  the  Bronze  Age  Zimmerman  (2009, 
2012)  argued  population  densities  in 
Central  Europe  between  0.6  and  1.8 
persons  per  100  km2 
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INFERRING  DEMOGRAPHIC  CHANGES 
FROM  THE  RADIOCARBON  RECORD 
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EVIDENCE  OF  AN 
INDIVIDUAL  ACTION 

(ISOTOPIC  EVENT) 
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DATASET  OF  14C  DATES 
(ARCHAEOLOGICAL  EVENT) 

■ > 

EVIDENCE  OF  A 
COLLECTIVE  ACTION 

Summing  a group  of  radiocarbon  dates  we  can  infer  demographic 
changes  in  a space-time  unit.  The  most  used  method  is  the  Summed 
Calibrated  Probability  Distribution  (SCPD). 
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POPULATION  STATIONARITY  BETWEEN  1800 
AND  800  BC  ON  A MACRO  SCALE 
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POPULATION  STATIONARITY  BETWEEN  1800 
AND  800  BC  ON  A MACRO  SCALE 


Since  we  do  not  want  to  count  twice  the  same 
demographic  signal  we  have  divided  between: 


Residential  areas 


Funerary  contexts 
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DIFFERENCES  IN  THE  TEMPORAL 
DISTRIBUTION  OF  COLLECTED  DATA 
BETWEEN  1800  AND  800  BC 
ON  A REGIONAL  SCALE 


Swiss  Plateau  Sum 
n=208 
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QjtGal  v4  2.3  Btonk  Ramsey  |2013),  f 5 


Padan  Plain  Sum 
n=109 


Calendar  date  (BC) 


Archaeological  sites  included  in  the  regional  SCPDs. 
The  Swiss  Plateau  (A),  the  Padan  Plain  (B), 
the  Massif  Central  (C)  and  the  French  Mediterranean 
coast  (D). 


GROWTH,  DIFFUSION  AND  THE 
ADOPTION  OF  INNOVATIONS 

ACROSS  TIME 
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SOCIAL  THEORY: 
CULTURAL  TRANSMISSION 


> The  Innovation 


Cremation  burials  and  new 
pottery  an  metallic  typologies 


Mediums  by  which  information  is 

^ Interaction  Channels  transmitted  to  or  within  a social 

system 


> Time 


> Space 


. , c . c "a  set  of  interrelated  units  that  are 

r The  Social  System  engaged  in  joint  problem-solving  to 

accomplish  a common  goal"  (Rogers 
2003) 


THEORETICAL  REMARKS 
The  S-shaped  curve 


13.5  % 34  % 34  % 


The  diffusion  of  innovation  according  to  Rogers  (2003). 
The  normal  distribution  is  in  blue  and  the  cumulative 
frequency  distribution  is  in  yellow. 
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Probability  density 


SCPDs  AND  HISTOGRAMS  OF  MEDIANS  OF  FUNERARY 
CONTEXTS  INCLUDED  IN  THE  EUBAR  DATABASE 

(we  did  not  take  into  account  dates  with  a standard  error  greater  that  100 

years) 


Inhumation  burials 


Cremation  burials 


OxCal  v4.2.3  Bronk  Ramsev  (2013):  r:5 

inhumation  Sum 
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OxCal  v4.2.3  Bronk  Ramsey  (2013):  r:5 


1800  1600  1400  1200  1000  800 

Calendar  date  (BC) 


ON  A REGIONAL 
SCALE 


SCPDs  of  14C-dated 
cremation  burials  from: 
the  North-East  of 
Iberian  Peninsula  (A), 
Southern  France  fB). 
Northern  Italy  (C)  and 
the  north  of  the  Alps 
region  (D). 
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5CPDS  AND  HISTOGRAMS  Or  MEDIANS  Or  14C-DAT6D 

CONTEXTS  WITH  FLUTED  POTTERY  INCLUDED  IN  THE 


EUBAR  DATABASE 


(we  did  not  take  into  account  dates  with  a standard  error  greater  that  100  years) 


JxCal  V4.2.3  Bronk  Ramsey  (2013);  r:5 

fluted  pottery  Sum 
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SCPDs  AND  HISTOGRAMS  OF  MEDIANS  OF  CONTEXTS  WITH 
FINDS  OF  HANDLE  WITH  VERTICAL  EXPANSION 


INCLUDED  IN  THE  EUBAR  DATABASE 

(we  did  not  take  into  account  dates  with  a standard  error  greater  that  100  years) 


OxCal  v4  2.3  Bronk  Ramsey  (2013):  r:5 
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NOT  ALL  THE  VARIABLES  CHARACTERIZE 

THE  NEW  PERIOD... 


OxCal  v4.2.3  Bronk  Ramsev  (2013):  r:5 
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DxCal  v4  2.3  Bronk  Ramsev  (2013):  r:5 
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Calendar  date  (BC) 


Probability  density  Probability  density 


SCPDs  AND  HISTOGRAMS  OF  MEDIANS  OF  RELIABLE 
CONTEXTS  WITH  DAGGERS  AND  METAL  KNIVES 


OxCal  v4  2.3  Bronk  Ramsey  (2013):  r:5 
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Calendar  date  (BC) 


Probability  density 


SCPDs  AND  HISTOGRAMS  OF  MEDIANS  OF 
FORTIFIED  SETTLEMENTS 
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FITTING  THE  EXPLANATORY  MODEL  TO 
ARCHAEOLOGICAL  DATA 


Generalized  logistic  distribution  fitted 
to  SCPDs  data 


y(t)  = a + 


K-A 


(1  + 


> cremation  burials  (A) 

> fluted  pottery  (B) 

> metal  knives  (C) 

> vases  with  handles 
with  vertical  expansion 
(D). 

The  green  dots 
represent  the  SCPDs, 
the  black  line  the  fitted 
curve  and  the  dashed 
blue  lines  the  95% 
confidence  interval. 

R2  is  the  coefficient  of 
detemination. 
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DATA  FROM  20  MACRO  REGIONS 


We  took  into  account 
20  regions  calculated 
using  Thiessen 

polygons. 

Each  region  constitutes 
a buffer  zone  around 
one  central  point. 


DATA  STRUCTURE:  Lower  and  Jpper  Boundaries 


QxCal  2 3 Bronk  Ramsey  I‘20l3x  r:S  InlCar  3 atnro&Dh&nc  curve  ■‘Re-imer  et  al  20131 
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Survival  time  (or  duration)  can  be  described  by  two  independent  parameters, 
a start  event  ta  and  a finish  event  tb. 
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ARCHAEOLOGICAL  DATA  FROM  20  GEOGRAPHICAL  REGIONS 


Region 

Funerary  ritual 

Lower  Boundary 

Upper  Boundary 

Ebro  Valley  Inhumation  1800  1519 


Ebro  Valley 

Cremation 

821 

800 

Barcelona 

Inhumation 

1800 

1334 

Barcelona 

Cremation 

997 

800 

W.  Pyrenees 

Inhumation 

1800 

1517 

W.  Pyrenees 

Cremation 

896 

800 

C.  Pyrenees 

Inhumation 

1800 

744 

C.  Pyrenees 

Cremation 

1358 

800 

E.  Pyrenees 

Inhumation 

1800 

890 

E.  Pyrenees 

Cremation 

960 

800 

Languedoc  Inhumation  1800  1435 
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W.  Pyrenees 
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1517 

W.  Pyrenees 
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800 

C.  Pyrenees 
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744 

C.  Pyrenees 
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1358 

800 

E.  Pyrenees 
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1600 

800 

Apennines 
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N.  Adriatic 
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N.  Adriatic 
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Vienna  Basin 

Inhumation 

1800 
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Vienna  Basin 
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CREATION  OF  TWO  SURVIVAL  TIME  TABLES: 


a)  for  measuring  the  survival  of  traditional  cultural  feature  (INHUMATION) 
b)  to  estimate  the  temporal  delay  in  adopting  the  cultural  change  (CREMATION) 

Tj . positive  number  indicating  the  number  of  calendar  years 


• SURVIVAL  TIME  TABLE  A: 
survival  of  BA  traditions  has  been 
calculated  as  the  difference  between 
time  0 (1800  BC)  and  the  median  of 
Bayesian  interval  for  the  last 
occurrences  of  inhumantions  in  a 
region  within  the  studied  period 


REGION 

INHUMATION 
SURVIVAL  TIME 
{numb  er  o f c alendar 
years  since  1 800  BC 
until  the  last 
occurrence) 

Ebro  Valley 

281 

Barcelona 

466 

W.  Pyrenees 

2 S3 

C.  Pyrenees 

1056 

E.  Pyrenees 

910 

Languedoc 

365 

Massif  Central 

651 

• SURVIVAL  TIME  TABLE  B: 

the  delay  in  adopting  the  funerary 
ritual  has  been  calculated  as  the 
difference  between  time  o (1800  BC) 
and  the  median  of  the  Bayesian 
interval  for  the  first  occurrence  of 
cremations  in  a region  within  the 
studied  period 


REGION 

CREMATION 
ADOPTION 
(numb  er  o f c alendar 
years  since  1S00 
BC  until  the  first 
occurrence) 

Ebro  Valley 

979 

Barcelona 

S03 

W.  Pyrenees 

904 

C.  Pyrenees 

442 

E.  Pyrenees 

S40 

Languedoc 

S15 

Massif  Central 

760 

THE  SURVIVAL  OF  INHUMATION  IN  BRONZE  AGE  EUROPE 


SURVIVAL  OF 

INHUMATION  BURIAL  AS  A 
TRADITIONAL  BURIAL 
PRACTICE 


t 


Y 


KAPLAN-MEIER  ESTIMATOR  OF  THE 


ANALYSIS  CONCERS  THE 
DURATION  OF  THE  LAST 
ARCHAEOLOGICAL 
EVIDENCES  OF 
INHUMATION 


SURVIVAL  FUNCTION 


Let  0 < t1  < t2  < t,  < ...  < tL  be  the  ordered  (continuous)  time  points  at  which  events  occur,  n,  is  the 
number  of  cases  at  the  risk  right  after  tl_1  and  dl  the  number  of  events  at  time  point  tv  Note  that  nl 
equals  nl_1  minus  the  number  of  events  and  the  number  of  censored  cases  in  time  interval  /.  Values  are 
plotted  using  a 95%  confidence  interval 
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THE  SURVIVAL  OF  INHUMATION  IN  BRONZE  AGE  EUROPE 


Survival  Plot 


From  1450  onwards,  the  probability  of  this  particular  tradition  surveying 
begins  to  decrease.  In  150  years,  probability  decrease  20%.  After  a 
relatively  short  period  without  further  changes  (from  1230  to  1120  BC), 
the  curve  record  a sudden  decrease  of  additional  20%.  Around  1120  BC 
inhumation  survived  in  60%  of  regions  from  the  Alps  to  the  Mediterranean. 
From  950  BC  onf  in  100  years  the  probability  of  inhumation  survival 
decreases  more  than  a 50%. 


THE  DELAY  OF  CREMATION  BURIALS  IN  BRONZE  AGE  EUROPE 


1800  BC 


Survival  Plot 


The  resulting  cumulative 
survival  curve  represents 
the  proportion  of  regions 
that  will  resist  to  cultural 
change  at  least  a given 
period  of  time  (delay). 


i K-M  Probability 
1 Lower  Limit 
: Upper  Limit 


Until  1700  BC  resistance  to  cultural  change  is  maximum  because 
inhumation  is  at  its  "floruit".  The  first  occurrence  of  cremation  can  be 
understood  as  a sequence  of  fast  but  relatively  unimportant  occurrences  of 
the  new  ritual.  At  1500  BCf  cultural  delay  was  high  enough  (80%).  The 
probability  of  cultural  change  resistance  maintains  for  another  100  years, 
after  which  decreases  quickly,  40%  in  200  years.  At  1200  BC  a majority  of 
regions  have  adopted  the  new  cultural  practice.  From  1000  BC  onwards,  the 
dynamics  of  cultural  change  increase,  so  that  at  750  BC  all  regions  from  the 
Alps  to  the  Mediterranean  have  adopted  cremation  rite. 


THE  RELATIONSHIP  BETWEEN  THE  SURVIVAL  OF  TRADITIONS  AND 
THE  ADOPTION  OF  NEW  CULTURAL  FEATURES 


Survival  Plot 


Log-Rank  test  using  Mantel  and 
Haenszel  statistic  produced  a p-value  of 
0.000,  which  suggests  that  there  are 
relevant  differences  between  both 
curves,  indicating  that  survival  of  old 
traditions  and  resistance  to  adopt 
innovations  were  not  temporally 
correlated  in  a linear  wav. 


K-M  Probability  Lower  Limit  Upper  Limit 


i Cremation  Cremation  Cremation 

i Inhumation  Inhumation  Inhumation 


The  probability  of  survival  of  the  tradition  was  always  higher  that  the  probability 
of  resistance  to  change.  Cremation  in  no  case  substituted  entirely  the  previous 
ritual.  The  moment  at  which  50%  of  studied  regions  maintained  the  old  tradition 
(1000  BC)  is  more  recent  than  the  moment  at  which  50%  of  regions  adopted  the 
new  funerary  ritual  (1200  BC).  Part  of  the  explanation  lies  in  the  quality  of  data 
used,  and  the  relevance  of  right-censoring  in  the  case  of  inhumation  burials. 
When  a majority  of  regions  (50%)  adopted  the  new  ritual  (1200  BC),  90%  of 
regions  still  used  the  traditional  practice. 


THE  PROBABILITIES  OF  CULTURAL  CHANGE  AT  THE  END  OF  BRONZE  AGE 


Hazard  function  of  adopting  a new 
ritual  at  the  end  of  Bronze  Age 


Use  of  a different  data  set 


time  0 as  the  median  of 
the  Bayesian  interval  for 
the  region  where 
cremation  first  occurred 


Region 

N.  Region 

Time  from  0 (Swiss  Plateau) 

Ebro  Valley 

1 

854 

Barcelona 

2 

678 

W.  Pyrenees 

3 

779 

C.  Pyrenees 

4 

317 

E.  Pyrenees 

5 

715 

La nguedoc 

6 

690 

Massif  Central 

7 

635 

Provence 

8 

159 

Lower  Rhone 

9 

475 

Upper  Rhone 

10 

352 

Alsace 

11 

422 

W.  Al  ps 

12 

717 

Swiss  Plateau 

13 

CO 

1675  BC=  0 223 


N-A  Cumulative  Hazard 


We  can  plot  the  cumulative  hazard  function  H(t),  or 
rather  an  estimate  of  it.  We  have  used  the  Nelson- 
Aalen  estimator: 


HNA{tt)  = E ^ 

3 H it-  L’ 

N-A  Cumulative  Hazard  Plot 


TIME 


Hasandi  cui'vssi 
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SPACE-TIME  DIFFUSION  OF 
THE  ADOPTED  INNOVATIONS 
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DETECTING  A SPACE-TIME  GRADIENT 


Expansive  phenomena  in  archaeology 

When  a system  expands  through  time,  we  can  foresee  a certain 
degree  of  dependence  between  locations,  and  this  dependence  is 
exactly  what  gives  an  appearance  of  unity  to  the  process. 


The  formal  conditions  for  an 
expansive  phenomenon  are 
the  existence  of  a spatial 
gradient  and  directivity,  which 
implies  a similarity  in 
neighboring  regions  (Tobler's 
Law). 
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MODELING  THE  FIRST  OCCURRENCE  OF 
CREMATION  BURIALS  BETWEEN  1800  AND  800 
BC  IN  PROTOHISTORIC  EUROPE 


Space-time  variations  in  the  first  occurence  of  cremation  burials  (medians 
of  the  calibrated  14C  dates).  Interpolation  has  been  calculated  using  a 
kriging  algoritm.  Countours  represent  differences  of  50  years.  Direccionality 
of  the  adoption  of  cremation  phenomenon. 

We  have  detected  an  average  rate  of  expansion  of  0.6-lkm/year. 


SPATIAL  AND  TEMPORAL  VARIATIONS  IN  THE  FIRST 
OCCURRENCE  14C  AND  TYPOLOGICALLY  DATED 

CREMATION  BURIALS 


Contours  represent  differences  of  50  years 
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DIRECTIVITY  OF  THE  2nd  MILLENNIUM  CREMATION 

ADOPTION  PHENOMENON 


The  map  contains  small  arrows  at  each  grid  node  pointing  down  the 
gradient,  that  is,  decrease  in  chronology:  from  places  where  cremations 
were  older  to  places  where  such  phenomenon  appears  to  be  more 
recent. 


THE  DIFFUSION  OF  VASES  WITH  HANDES  WITH 
VERTICAL  EXPANSION  IN  THE  MIDDLE  AND  LATE 

BRONZE  AGE 
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THE  DIFFUSION  OF  VASES  WITH  HANDES  WITH 
VERTICAL  EXPANSION  IN  THE  MIDDLE  AND 
LATE  BRONZE  AGE:  FROM  EAST  TO  WEST 
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Map  showing  the  directivity  of  the  vases  with  handles  with 
vertical  expansion  adoption  phenomenon 


NOT  ALL  THE  DIFFUSION  PROCESSES  FOLLOW  THE 

WAVE  OF  ADVANCE  MODEL... 


Spatial  and  temporal  variations  in  the  first  occurrence  of  contexts 
characterized  by  the  presence  of  fluted  pottery  included  in  the 
EUBAR  database  (medians  of  the  calibrated  radiocarbon  dates). 
Contours  represent  differences  of  50  years. 
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CONCLUSIONS: 

Cultural  Change  at  the  End  of  Bronze  Age 
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Cal  Pal  (Version  March  2007) 


452  radiocarbon  dates,  covering  a span  of  ~ 3000  years,  between 
9700  and  6700  cal  BC  (From  the  PPNA  to  the  Early  Pottery 
Neolithic  [PN]),  shows  a significant  reduction  in  the  volume  of 
the  archaeological  signal  around  8000  cal  BC,  indicating  that 
the  apparent  uninterrupted  sequence  of  settlements  from  PPNA 
to  PN  cannot  be  taken  as  unequivocal  evidence  of  settlement 
continuity 


1 

DEPARTAMENTS  de  PREHISTORIA 
i CIENCIES  DE  L’ANTIGUITAT  ^ v 


UflB  0 

IJniviTsitat  Autonoma 

dc  Barcelona  wtD«Mamu!^u 


Proximo  oriente 


CalPal  (Version  March  2007) 


l-yS 

Jerfel  Ahmar 

95%  Peak  9680  cal  BC  - 8430  cal  BC  Span  =1250 
68%  Peak  9430  cal  BC  • 8790  cal  BC  Span  = 640 
50%  Peak  9350  cal  BC  • 8890  cal  BC  Span  = 460 

N = 24  A ™J\ 

fljV1  | 

Mureybet 

95%  Peak  9720  cal  BC  • 7960  cal  BC  Span  =1760 
68%  Peak  9400  cal  BC  - 8450  cal  BC  Span  = 950 
50%  Peak  9290  cal  BC  - 8590  cal  BC  Span  = 700 

: N = 20 

y ■ n ■■  i ii  m i V 

Qaramel 

95%  PeaklOOOl  cal  BC  - 8690  cal  BC  Span  =1320 
68%  Peak  9720  cal  BC  • 8950  cal  BC  Span  = 770 
50%  Peak  9620  cal  BC  - 9100  cal  BC  Span  = 520 

: N=T 

u 

'Abr  3 

95%  Peak  9290  cal  BC  • 8860  cal  BC  Span  = 430 
68%  Peak  9250  cal  BC  - 8970  cal  BC  Span  = 280 
50%  Peak  9240  cal  BC  - 9100  cal  BC  Span  = 140 

:"s  y 

Djade 

95%  Peak  9180  cal  BC  - 8160  cal  BC  Span  =1020 
68%  Peak  8890  cal  BC  ■ 8340  cal  BC  Span  = 550 
50%  Peak  8740  cal  BC  - 8380  cal  BC  Span  = 360 

Gobekli 

95%  Peak  9670  cal  BC  • 8400  cal  BC  Span  =1 270 
68%  Peak  9540  cal  BC  • 8790  cal  BC  Span  = 750 
50%  Peak  9460  cal  BC  • 8890  cal  BC  Span  = 570 

:n-”  "a| 

J , p 1 1 1 

u no — 

Nevali  £ori 

95%  Peak  8600  cal  BC  - 7200  cal  BC  Span  =1400 
68%  Peak  8490  cal  BC  - 7800  cal  BC  Span  = 690 
50%  Peak  8440  cal  BC  - 7940  cal  8C  Span  = 500 

'!& 

Ain  el  Kerkh 

95%  Peak  8650  cal  BC  - 5820  cal  BC  Span  =2830 
68%  Peak  8480  cal  BC  - 6440  cal  BC  Span  =2040 
50%  Peak  8420  cal  BC  - 6630  cal  BC  Span  =1790 

N = 26 

PL  "1 

Abu  Hureyra 

95%  Peak  8080  cal  BC  - 6740  cal  BC  Span  =1340 
68%  Peak  7670  cal  BC  - 7090  cal  BC  Span  = 580 
50%  Peak  7560  cal  BC  - 7170  cal  BC  Span  = 390 

: N = 19° 

/ 

- 

llll  Ill  II  C~r-7||VT\ 

Halula 

95%  Peak  81 10  cal  BC  - 6310  cal  BC  Span  =1800 
68%  Peak  7810  cal  BC  - 6820  cal  BC  Span  = 990 
50%  Peak  7720  cal  BC  • 7210  cal  BC  Span  = 510 

: N = 26 

mini  mi  iW-f-'Tv 

Akargay 

95%  Peak  7790  cal  BC  - 6190  cal  BC  Span  =1600 
68%  Peak  7430  cal  BC  • 6790  cal  BC  Span  = 640 
50%  Peak  7360  cal  BC  - 6940  cal  BC  Span  = 420 

: N = 11 

LU 111  III  V 

Gritille 

95%  Peak  8200  cal  BC  - 6660  cal  BC  Span  =1540 
68%  Peak  7910  cal  BC  • 6910  cal  BC  Span  =1000 
50%  Peak  7780  cal  BC  - 7030  cal  BC  Span  = 750 

: n=T 

Li 

to  ~ 1 

Hayaz/Kumartepe 

95%  Peak  7490  cal  BC  - 6640  cal  BC  Span  = 850 
68%  Peak  7410  cal  BC  - 6750  cal  BC  Span  = 660 
50%  Peak  7350  cal  BC  - 6820  cal  BC  Span  = 530 

N = 174 

Sabi  Abyad 

95%  Peak  71 10  cal  BC  - 5890  cal  BC  Span  =1220 
68%  Peak  6820  cal  BC  • 6130  cal  BC  Span  = 690 
50%  Peak  6630  cal  BC  - 6220  cal  BC  Span  =410 

: N=T 

i 

Damishliyah 

95%  Peak  7050  cal  BC  ■ 6430  cal  BC  Span  = 620 
68%  Peak  6910  cal  BC  ■ 6490  cal  BC  Span  = 420 
50%  Peak  6820  cal  BC  • 6520  cal  BC  Span  = 300 

: N = 25 

J 

Oh"  1 

Bouqras 

95%  Peak  7500  cal  BC  • 6290  cal  BC  Span  =1210 
68%  Peak  7380  cal  BC  - 6680  cal  BC  Span  = 700 
50%  Peak  7280  cal  BC  - 6750  cal  BC  Span  = 530 

' n=T 

A 

a '1 

Sinn 

95%  Peak  7750  cal  BC  - 6930  cal  BC  Span  = 820 
68%  Peak  7640  cal  BC  - 7160  cal  BC  Span  = 480 
50%  Peak  7480  cal  BC  ■ 7210  cal  BC  Span  = 270 

: N = 10 

/L 

1 i 

Shir 

95%  Peak  7040  cal  BC  • 6450  cal  BC  Span  = 590 
68%  Peak  6920  cal  BC  • 6510  cal  BC  Span  =410 
50%  Peak  6830  cal  BC  - 6550  cal  BC  Span  = 280 

' N = 10 

, A.  : 

Seker  al  Aheimar 

95%  Peak  7070  cal  BC  - 6400  cal  BC  Span  = 670 
68%  Peak  6940  cal  BC  ■ 6530  cal  BC  Span  =410 
50%  Peak  6850  cal  BC  - 6580  cal  BC  Span  = 270 

: NT 

JU 

Mezraa  Teleilat 

95%  Peak  7060  cal  BC  - 6520  cal  BC  Span  = 540 
68%  Peak  7000  cal  BC  - 6670  cal  BC  Span  = 330 
50%  Peak  6960  cal  BC  - 6730  cal  BC  Span  = 230 

10000  9000  8000  7000  6000  [calBC] 


U 


Laboratori  d’Arqueotogia  Quantitativa 


Tw&fihi  p 

PTiilippon 

LC2i  100  j 

i **  - 

J5-  s o-l 


[%I  -SHI 
-10 
im 


\jm 


Rg  Id.  Palaeocl  innate  p noxies  fro  m both  the  Eastern  Mediterranean  and  the  North  Atlantic.  Note  Ihe 
syrKhtforriy  telvwee n 1)  Ire  1G.£  cal  Hy  BP  RCG  amd  re  mailed  decrease  n re  aichaedlogicail  signal  n re 
nerfieni  Levanl , and  2\  line  boosl  observed  n Ihe  aic  naeotgraS  signal  and  ine  beginning  o\  re  ealty 
HolQQene  sapropei  Si  period  (modified  kom  rel . |7Qfl. 


What  does  it  mean 
Spatio-Temporal  ANALYSIS? 


Distinction  and  separation  of  parts  of  spatio-temporal 
ordering  of  archaeological  observables  in  order  to  know 
their  relational  structure, 

Qualitative  and  quantitative  examination  of  spatio- 
temporal  relationships 

Solving  the  where,  when  and  why  problems  in 
archaeology. 


